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Operators on unbounded domains may acquire eigenvalues that are embedded in the essential spectrum. 
Determining the fate of these embedded eigenvalues under small perturbations of the underlying operator is 
a challenging task, and the persistence properties of such eigenvalues is linked intimately to the multiplicity 
of the essential spectrum. In this paper, we consider the planar bilaplacian with potential and show that the 
set of potentials for which an embedded eigenvalue persists is locally an infinite-dimensional manifold with 
infinite codimension in an appropriate space of potentials. 

> 

On 

1 Introduction 

O 



Determining the dependence of the spectrum of operators on perturbations is an important issue that is of 



o 

relevance in many applications. Of course, much is known in this direction: the persistence of point eigenvalues 
and the behavior of the essential spectrum under small bounded perturbations, for instance, have been analyzed 
comprehensively, and we refer to [10] for many results along these lines. Here, we consider differential operators 
that are posed on unbounded domains and are interested in the interaction between eigenvalues, with proper 
cigenfunctions in the underlying domain of the operator, and the essential spectrum. More precisely, we study 
the fate of eigenvalues that are embedded in the essential spectrum under small perturbations of the operator. 
Typically, such eigenvalues will disappear under generic perturbations of the potential, and it is therefore of 
interest to determine the class of perturbations for which an embedded eigenvalue persists. For the bilaplacian 
on cylindrical domains, we showed in our previous work [7] that the set of perturbations for which an embedded 
eigenvalue persists is an infinite-dimensional manifold of finite codimension. Furthermore, we showed that 
the codimension of this set is given by the multiplicity of the essential spectrum, defined as the number of 
independent continuum eigenfunctions or, more rigorously, via the spectral resolution of the Fourier transform 
of the bilaplacian (see e.g. [2, Definition 2 in §85]). In this paper, we continue the investigation that we began 
in [7] and consider the bilaplacian posed on the plane: the challenge is that the essential spectrum of the planar 
bilaplacian has infinite multiplicity. Thus, wc may expect that the set of potentials for which an embedded 
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eigenvalue persists is an infinite-dimensional manifold of infinite codimcnsion, and this is indeed what we shall 
prove for an appropriate class of potentials. 

Before stating our results, we briefly outline why embedded eigenvalues are of interest. Our first motivation comes 
from quantum mechanics: the eigenfunctions associated with eigenvalues of an energy operator correspond to 
bound states that can be attained by the physical system modelled by the energy operator. If such an eigenvalue 
is embedded in the essential spectrum, then its fate under perturbations of the potential determines whether the 
associated bound states persists or not (see [9, 15] for examples). A second example comes from inverse scattering 
theory, where eigenvalues correspond to coherent soliton structures of the underlying intcgrable system, while 
the essential spectrum describes radiative scattering behavior. Thus, bifurcations of solitons arc reflected by 
the disappearance or persistence of embedded eigenvalues [12, 13]. Finally, embedded eigenvalues provide a 
common mechanism for the destabilization of travelling waves in near-intcgrablc Hamiltonian partial differential 
equations, and we refer to [16] for further background information and pointers to the literature. 

As mentioned above, we focus in this paper on the persistence of embedded eigenvalues for the planar bilapla- 
cian. Our primary reason for considering the bilaplacian is that this operator is complex enough to exhibit the 
underlying difficulties, while not adding technical complications that have nothing to do with the issue we are 
interested in. In other words, the planar bilaplacian provides a useful paradigm for the issues that we expect to 
encounter for other more complicated differential operators. Note also that the applications we mentioned above 
all involve selfadjoint operators, a feature shared by the bilaplacian. 

We now describe the precise setting that we consider. Let ro > 0, and assume that 9 £ C§°(B ro (0); R) is a radially 
symmetric potential. Hence, we use polar coordinates (r, <p), write 9 = 6(r), and consider the multiplication 
operator on L 2 (R 2 ) (also denoted by 9) defined by 

[6u](r,ip) := 9(r)u(r,ip). 

We define C := A 2 + 9 on L 2 (R 2 ), where A 2 is the bilaplace operator which is densely defined on L 2 (M. 2 ) with 
domain fJ 4 (R 2 ). It is known that the spectrum of A 2 is cr(A 2 ) = [0, oo). Since 9 has compact support, the 
essential spectra of C and A 2 coincide, and so a c (C) — [0, oo). We assume that 9 is chosen so that C has a simple 
positive eigenvalue Ao: 

(Al) £ has an eigenvalue Ao > of multiplicity 1. 

We are mainly interested in the case where Ao is an embedded eigenvalue, i.e. when Ao>0, since when Ao is 
isolated from the rest of the spectrum, the persistence of eigenvalues is well known, [10, pp. 213-215]. We also 
exclude the case Ao — which lies on the boundary between spectrum and resolvent set. We denote by u*(r, <p) 
the eigenfunction associated with the embedded eigenvalue Ao- Since 9 is radially symmetric and the laplacian 
A is invariant under rotations of the underlying cartesian coordinates, we see that the functions u*(r,(p + <po) 
are, for each fixed tpo, also eigenfunctions of L belonging to the eigenvalue Ao- The simplicity of Ao required in 
assumption (Al) therefore implies that is a radial function, and we henceforth write it* = it*(r). It is clear 
by existence and uniqueness of solutions of ODEs that (r) cannot vanish for all r > ri, and so we assume that 

(A2) n > ro is such that u*(ri) 7^ 0. 

Lemma 1 below shows that our hypotheses can be satisfied. We now perturb the potential 9 by potentials p in 
the weighted L 2 -spacc TZ := L 2 ([0, n], i? 1 / 2 (S' 1 ), r dr) of functions that map the interval [0,ri] into H 1 ^ 2 (S 1 ), 
where the interval [0,ri] is the domain of the radial variable r, while H 1 / 2 (S 1 ) describes the dependence on the 
angular variable <p. Our main result is as follows. 

Theorem 1. Let < ro < ri, 6 € Cff (B r<1 (0); R) be radially symmetric, and assume that (Al) and (A2) hold. 
Then there exists 6 > and a neighbourhood O of in TZ = L 2 ([0, n], if 1/2 (S 1 ), r dr) such that the set 

TZemb '■= {p S O; C + p has an embedded eigenvalue in (Ao — 5, Ao + 5)} 

is a smooth manifold in TZ of infinite dimension and codimension. 
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Before commenting on the ideas behind the proof of Theorem 1, we illustrate that our hypotheses can be met. 

Lemma 1. There exists a smooth radial potential 6(r) with compact support such that C = A 2 + 9 satisfies 
Hypothesis (Al). 

Proof. Let K${r) denote the modified Bessel function of the second kind and define a smooth, strictly positive 
function Uo(r) via 

1 < r < 1 

K (r) 2<r 

together with a smooth interpolation in the intermediate region r £ [1,2]. Note that uq decays to exponentially 
as r — v oc and can be chosen so that Uo{r) > for all r. Thus, the radial potential 



uo(r) 



-I(-A 2 + lH = -(A + l)(-A + l) Uo = 1 l~ r ~ l 
u u 2 < r 



is well-defined, smooth, and has support in [0, 2] since Kq{t) satisfies (—A + l)A'o = 0. Furthermore, we have 

Cu = (A 2 + 6)u Q = A 2 u + — [(-A 2 + l)uo\u = u , 

and uo is a positive radial eigenfunction belonging to the embedded eigenvalue Ao = 1 of C. 

It remains to show that Ao = 1 is simple. Using the radial symmetry of 9, the results presented in the rest of 
this paper imply that it suffices to show that the equation 

2 



1 k 2 

dl + - r d r - ^ 



u={l-ff)u for r£ (0,2) (1) 
u r (0) = u rrr (0) = 0, u{r) = K k {r) for r > 2 



does not have a solution u(r) for each integer k > 1. We will now outline why (1) will not have solutions for k > 1 
provided 9 is modified appropriately but omit the straightforward details. Using variations of parameters, it can 
be shown that (1) cannot have solutions for k 3> 1. If it does have solutions for some or all of the remaining 
finitely many integers k > 1, then we can modify the potential 9 to remove these solutions while retaining the 
eigenfunction for k — 0. Indeed, any solution of (1) for k > 1 is of the explicit form u(r) = r 1, for some integer 
I = £(k) > 1 since 6(r) — 1 = for < r < 1. Replacing uo in the above construction of 9 by uo + evo for 
bounded functions vq with support in (1, 1) and using the necessary expressions (64) derived in section 7 for the 
persistence of eigenvalues, it is then not difficult to see that any nonzero choice of i>o > removes the solutions 
of (1) for k > 1 for 0«e«l. □ 

The idea for proving Theorem 1 is to characterize embedded eigenvalues as roots of a regular function, since 
such a characterization would allow us to use the implicit function theorem. As it appears difficult to find a 
functional-analytic characterization of embedded eigenvalues, we pursue here a dynamical-systems formulation 
similar to that used in our precursor work [7] for the bilaplacian on cylinders. The eigenvalue problem can be 
written as a system of differential equations in the radial evolution variable r posed on an appropriate function 
space X of functions that are defined in the angular variable tp. The issue is that this system is ill-posed in the 
sense that, for given initial data, solutions may not exist. Using a similar approach as in Scheel [18], we will 
show, however, that this dynamical system has an exponential dichotomy: there are two infinite-dimensional 
subspaccs X cu and X s of X at r = n so that solutions with initial data in X cu exist and stay bounded for 
r < ri, while solutions with data in X s exist and decay as r — > oo. The intersection of these spaces corresponds 
to eigenfunctions of the underlying operator, and our goal is therefore to characterize those perturbations for 
which this intersection is nontrivial. We show that there are infinitely many conditions that characterize such 
intersections and prove that we can solve them using an implicit function theorem. A key issue is the space for 
the perturbation p. For the conditions of the implicit function theorem to be satisfied, the space for p needs to be 
i 2 ([0, r{\] H 1 ^ 2 (S 1 ), r dr), a space with very low regularity. This low regularity forces us to work with different 
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function spaces for r < r± (where p has its support) and for r > r\ (where we have an explicit formulation of 
the solutions of the system in terms of Bessel functions) , and so we need to take extra care when matching the 
solutions at r = r\ . 

The rest of this paper is organized as follows. In section 2, we introduce the spatial-dynamics formulations of the 
eigenvalue problem. In sections 3 and 4, we prove the existence of exponential dichotomies for the bilaplacian 
and for the operator C, respectively, near the core r = 0. We then construct dichotomies for C in the far field 
for r > 1 in section 5 and discuss similar properties for the adjoint spatial dynamical system in section 6. 
These results are then used in section 7 where we match the solutions from the core and the far field by using 
Lyapunov-Schmidt reduction and prove Theorem 1. The paper is concluded with suggestions for extensions and 
some open problems. 



2 Spatial-dynamics formulation 

If A is an eigenvalue of C + p, then there exists u £ H 4 (M. 2 ) such that 

A 2 u + (6 + p)u = Xu. 
Let r3 > T2 > max(l,r"i). We introduce a new radial variable 



s{r) 



logr if r < r-2, 
if r > r 3 , 



(2) 



(3) 



and for r £ (r 2 , r^), we define s such that s £ C°°(IR + ; K) is strictly increasing. Note that this implies that there 
exist constants c and C such that < c < C and c < s'(r) < C for every r 2 < r < r 3 . We define 9 and p by 
9{s{r)) = 9(r), etc. Since s is an increasing function, it is invcrtible, and we denote the inverse function by r(s). 
Let Sj := s(rj), j = 1, . . . , 3. Under the coordinate transformation (3), the space 1Z transforms into the space 1Z 
given by 

£:=L 2 ((-oo,Si]; H 1 ? 2 ^ 1 ), e 2s ds), 
that is, the weighted L 2 space with values in iJ 1/ ' 2 (S' 1 ) and weight e 2s . 
Setting v = Am, equation (2) is equivalent to the system 



An = v, 
Av = (A - 

where in the variables s and ip, the Laplacian is given by 

1 



p)u, 



(4) 



A = 



r'(s) 2 



/r'(s) r"(s)\ d 
ds 2 \ r(s) r'(s) J ds 



is) 



d 2 

d^ 2 



Rewriting this intermediate system as a first order system, with u\ 
denotes differentiation with respect to s, we obtain a system of the form 

U'(s) = A(s;X,p)U, 



u, U2 — u', U3 = v and 114 = v', where ' 

(5) 



where A(s; A, p) is given by 



A(s;X,p) 





r'(s) 2 q2 
■ r( S y ° 





r ' (s) r (a) 
r f (s) r(s) 







r'{s) 2 




\[X-~e-p)r'{sf 







(6) 



r '( s ) 2 g2 r"(s) _ r'(s) . 
r(s) 2 r'(s) r(s) / 



where d denotes differentiation with respect to p, i.e., d 
for s < S2 or s > S3; see sections 3 and 5. 



The expression for A(s; A, p) simplifies significantly 
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The perturbation p G 1Z is, in general, not continuous or even bounded, so we need to study weak solutions of 
(5). Let 

X = H 2 ^ 1 ) x H 1 (S 1 ) x H 1 (S 1 ) x L 2 (S 1 ); 
Y = H^S 1 ) x H^S 1 ) x H^S 1 ) x H 1 (S 1 ). 

Definition 1. Let J be an interval o/R. A function U : J — > X is a weak solution of (5) in J if 

1. UeLl c (J;Y)nHl oc (J;X), 

2. for every V G C§°(J; X) we have 

- f U{s)V / (s)ds = [ A{s;X,p)U(s)V{s)ds. 
J. j J j 

Lemma 2. Let A 6 R. The eigenvalue equation (2) has a solution u G _ff ; 4 oc (R 2 ) if and only if (5) has a weak 

solution u g ff/ oc (K; x) n £ 2 oc (m ; y ) n £°°(M-;^Q. 

Proof. Suppose that u G Hf oc (R 2 ) is a solution of (2), and let U := (u, u', Au, (Aw)') t , where ' denotes differen- 
tiation with respect to s. 

We first consider u as a function of r, and let -Br(O) be a ball centered at 0, with R any positive radius. Then 

ue H 1 ((0,R);H 3 (S 1 ),rdr)nL 2 ((Q,R);H i (S 1 ),rdr) 

C H 1 {(0,R);H 2 (S 1 ),rdr) n L 2 ((0, i?); ff 3 ^ 1 ), r dr), 

v! = G ffVfO, JZ); ffW), r dr) n £ 2 ((0, iZ); ff 3 ^ 1 ), r dr) 
dr 

C i/^O, fi); i?^ 1 ), r dr) n L 2 ((0, i?); iJ^S 1 ), r dr), 
Au G ^((0, iZ); iJ^S 1 ), r dr) n L 2 ((0, iZ); H 2 (S 1 ),rdr), 

(Au)' = r '^4^. g tfVfO, iZ); ^(S 1 ), r dr) n L 2 ((0, R);H 1 (S 1 ), r dr), 
dr 

where r' = dr/ds — r for r < r 2 . By the Sobolev embedding theorem, {/ G C([0, R]; X), and so U(s(r)) has a 
limit as r — > 0+, or equivalently, as s — » — oo. Hence, viewing {/ as a function of s, U G L°°(R_;X). We also 
see that C/ G i?j oc (R; X) n L 2 oc (R; F). It is clear from the construction that U is a weak solution of (5). 

Conversely, let U G ff/ oc (R; X) n Lf oc (R; Y) n L°°(]R_; X) be a weak solution of (5), and let u = V X - Viewing u 
as a function of r rather than of s, it is clear that u G i?^ c (R+; i? 2 (5' 1 )) C C((0, oo); C(S' 1 )), and so by (4), 

A 2 u = (A - 9 - P )u G L 2 oc (R+; l^S 1 ), r dr) C L 2 oc (R + ; L 2 {S 1 ), r dr) = L 2 oc (M 2 \ {0}). 

Since we also have 

ue L oo ([0,r 2 };H 2 (S 1 )) C L°°([0, r 2 ]; ^(S 1 )), 
A-ff-pGL 2 ([0,r a ];ir 1 / 2 (5 1 ),rdr) C L 2 ([0, r 2 ]; L 2 ^ 1 ), r dr), 
we see that A 2 it = (\-6-p)u G i 2 ([0,r 2 ];i 2 (5 1 ),rdr) = L 2 (B ? . 2 (0)). We have proved that it G H 4 (B r2 {0))n 

Since it is also clear that u solves (2), the proof is complete. □ 

Note that a weak solution satisfies U G C(M.;X) (see e.g. [8, p. 286]), and so the following definition for an 
exponential dichotomy makes sense (see also [6] for the standard definition for ODEs and [14] for an extension 
to PDEs): 

Definition 2. Let J be an unbounded subinterval o/R. We say that equation (5) has an exponential dichotomy 
in X on J if there exists a family of projections P(s) for s G J such that for any s G J, P{s) G C(X), 
P(s) 2 = P(s) and P(-)U G C(J;X) for every U G X, and there exist constants K > and k s < k u with the 
following properties: 
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(i) For each t G J and t/el there exists a unique weak solution Q s (s,t)U of (5) defined for s > t, s,t G J 
such that <$> s (t,t)U = P(t)U and 

\\$ s (s,t)U\\ x < Ke K3 ( s -V\\U\\ x 

for all s > t, s,t G J . 

(ii) For each t G J and i/iere exists a unique weak solution & u (s,t)U of (5) defined for s < t, s,t G J 
sucft ttat <S> u {t,t)U = (I- P(t))U and 

\\$ u {s,t)U\\x < Ke K ^ s -^\\U\\ x 

for all s < t, s,t G J . 

(Hi) The solutions $ s (s,t)U and & u (s,t)U satisfy 

$ s (s,t)U G RanP(s) for every s > t, s,t G J, 

<I> u (s,t)J7 G kcr P(s) for every s < t, s,t G J. 

We also need the definition of time-dependent exponential dichotomy, which will be used for J = [si,oo) and 
with X s := H 1 x L 2 x H 1 x L 2 with the s-dependent norm 

\\U\\ 2 X s := -^IKHffi(si) + \\ u i\\h(si) + IM|i=(,si) + ^I|w3||ffi (s i) + IMll^si) + H^IU^si), 
where itj are the components of U, j = 1, . . . ,4. 

Definition 3. Lei J 6e an unbounded subinterval of K. FKe say i/iai equation (5) /ias a time- dependent expo- 
nential dichotomy in X" on J if there exists a family of projections P{s) for s G J suc/i £/ia£ /or any s G J, 
.P(s) € £(A" S ), -P(s) 2 = P(s) and P(-)U G C(J;Af) /or every [/ G X, and there exist constants K > and 
k s < k u with the following properties: 

(i) For each t G J and U G X there exists a unique solution (fr s (s,t)U of (5) defined for s > t, s,t G J smc/i 
tfiat $ s (t,t)U = P(t)U and 

\\$ s (s,t)U\\ x « < Ke KS ( s -V\\U\\ xt 

for all s > t, s, t G J. 

(mJ For eac/i i G J and U G X there exists a unique solution <& u (s,t)U of (5) defined for s <t, s,t G J suc/l 
tfiai $ u (t,t)U = (I- P(t))U and 

\\$ u (s,t)U\\ X s < Ke« u ( s -V\\U\\ xt 

for all s < t, s,t G J . 

fmj T/ie solutions $ s (s,t)U and $ u (s,t)U satisfy 

& s (s,t)U G RanP(s) /or every s > t, s,t G J, 

& u (s,t)U G kcrP(s) /or every s < t, s,t G J. 

In the following sections, we will consider the intervals 

J_ := (— oo, Si] and J + := [si,oo), 

and show that the system (5) has an exponential dichotomy on J_ and a time-dependent exponential dichotomy 
on J + . 
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3 Dichotomies for the system at — oo 



For s < si, r(s) = e s , and hence the system (5) is given by 

u[ = u 2 , 

u' 2 = -d 2 ui + e 2s u 3 , 
u' 3 = Ua, 



(7) 



u' 4 = (A - 9(a) - p(s, •)) e 2s Ul - d 2 u 3 . 
In the limit as s — > — oo we have the system 







f 


1 





»\ 




M 


u' 2 




-d 2 




























1 






w 




^ o 





-d 2 






w 



or 



U'=A-U. 



(8) 



We expand C/ = (mi, U2, U3, U4) 7, as a Fourier series in the <p variable, and denote the fcth Fourier coefficient by 
Uk(s). For j <G R, we define the weighted / 2 spaces I 2 with norm defined by 



||{a fc } feez |||:=£(l + fc 2 )' 
fcez 



Ofc 



The function space induced by X is 



X :=l\ x x l\ x Z 2 . 
The system (8) decouples in the Fourier space and for £ Z we have 

^(s)=l_(fc)C/ fe ( S ), 

where 



(9) 
(10) 



A.{k) := 



^0 10 0^ 
k 2 
1 
\o fc 2 0/ 

The eigenvalues of A_ (fc) are ±|fc|, and for fc 7^ both eigenvalues have geometric multiplicity 2. The eigenvectors 
for k 7^ are (±l/fc 2 , l/|fc|, 0, 0) T and (0, 0, ±l/|fc|, 1) (we normalize the eigenvectors so that their X norm is 
approximately constant and bounded away from as fc — > 00) . Let 

1/fc 2 ^ 
l/|fc| 

■1/1*1 i/l^l 

1 / 

and 



D, := 



M k := 



/-1/fc 2 

1/1*1 







1 





f-\k\ 





\ 





-|fc| 











1*1 


V 





\k\) 



so that A-(k) = M k D k M^ x for fc ^ 0. Note that 



T 1 





/-fc 2 


1* 


1 








2 


fc 2 


1* 




I 






0\ 

-1*1 1 



1*1 1/ 
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For k = 0, the eigenvalue has algebraic multiplicity 4 and geometric multiplicity 2, so A_(0) is not diagonaliz- 
able. Note however that ^4-(0) is already in Jordan normal form. We define Mq = I and Do = A_(0). 

Lemma 3. The operator A_ : X — > X is a closed densely defined operator with spectrum a(A^) = Z. 

Proof. Recall that X = H 2 ^ 1 ) x H^S 1 ) x H^S 1 ) x I?{S l ), and Y = H 3 ^ 1 ) x i? 2 ^ 1 ) x i? 2 ^ 1 ) x H^S 1 ). 
It is easy to check that the domain of A- is Y, which is dense in X. To see that A_ is closed, let Uj £ Y be 
such that Uj — > U in X and A—Uj — > / in X. We write Uj = (ui 1 j,U2,j,U3 i j,Ui l j) T etc. By the definition of A_ 
wc have 



Ul,j 
U2,j 
1*3,3 
Ui,j 



— > U\ 



U3 
"4 



in H 2 , 
in H 1 , 
in if 1 , 
in L\ 



while 



"2,3 
U 4 ,j 



fi 

Pi 
h 
h 



in H 2 , 
in i? 1 , 
in H 1 , 
in L 2 . 



if 2 D If 3 , and since limits (in 



It follows that U2 = fi £ H , and that wi.j converges in H . Since u\ t j Ui m 
H 2 ) are unique if they exist, we also have u\.j — > U\ in H 3 , and so —d 2 u\ — fii- It follows in exactly the same 
way that m = fs £ H 1 , that us £ H 2 and that —d 2 u 3 = f±. This shows that U E Y and A-U = F, and so 
A- : X -> X is closed. 

The operator A^ : X — > X induces an operator A_ : X — > X defined by 

(A_U) k :=A_(k)U k . 
Then A_ is a densely defined operator on X with domain Y := I 2 x Z| x Z| x Z 2 . 

It is clear that (A_ — fil) : X — > X has a bounded inverse if and only if (A- — fil) : X — > X has a bounded 
inverse. It is also clear that fc <G cr(^4_) for k £ Z. To prove that there are no other points in the spectrum of 
A_ , let /i G C \ Z. 

Define M : I 2 x Z 2 x I 2 x Z 2 -> Z 2 x I 2 x Z 2 x I 2 by 

(Mft) fc = Mfc&fc, 

and note that M is a linear homeomorphism between these spaces. Define also the unbounded operator D on 
l 2 xl 2 xl 2 x I 2 by 

(DU) k = D k U k . 

Note that D is a closed densely defined operator with domain l\ x l\ x Z 2 x Z 2 , and that cr(£>) = Z. 
If /x e C \ Z, then 

(1_ - fiiy 1 = M(D - liiy^M' 1 
It is now easy to see that (A- — : X — » X is bounded, and consequently also (A_ — : X — » X. □ 

Having established that the spectrum of A^ consists exactly of its eigenvalues, we define the (generalized) 
spectral projections P s , P c , P u in X, corresponding to the negative, the zero and the positive eigenvalues of A-, 
respectively. Let X s = P S X, etc. so that X — X s Q)X C Q)X U , where X s and X u arc infinite-dimensional whereas 
X° is four-dimensional. We also define corresponding spectral projections P%, P£ of A-(k), in the spaces X k , 
k £ Z \ {0} and note that if U = Y. k & U k e ih £ X, then 

P£U k e ik \ 



P S U 



r k 

fcez\{o} 



fcez\{o} 
P C U = U . 
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Lemma 4. The operator A- possesses an exponential dichotomy in X on J_ = (—00, s±] with constant K and 
rates k s = and k u = 1, and another exponential dichotomy in X on J_ with constant K and rates k s = — 1 
and k u = 0. 

Proof. Let r\ £ (0, 1) be arbitrary. We apply Lemma 2.1 of [14] for the operators — nl and A_ + nl, and 
obtain exponential dichotomies with constant K and rates n s = —r\ and k" — 1 — and k s = — 1 + 77 and k u = 77, 
respectively. The existence of exponential dichotomies for A_ with rates k s ~ — 1 and k u = 0, and k s = and 
K 11 = 1, respectively then follows by using the transformation V = e ±v 'U. 

We only consider the operator A_ — r/I, since the proof for A- + rjl is similar. The result follows from Lemma 
2.1 of [14] if we can verify condition (HI) of [14] for the operator A^ — namely 



(HI) Suppose that there exists a constant C > such that 

\\(A- -ifiiy 1 ] 

for every pet. 



C(X) 



< 



C 



i + |m| 



As in the proof of Lemma 3, it suffices to prove that there exists a constant C such that 



D — r/I — i\xl 



U 



< 



C 



l 2 xl 2 xl 2 xl 2 



1 + |m| 



1^1 



l 2 xl 2 xl 2 xl 2 



for every U £ I 2 x I 2 x I 2 x I 2 . Note that for fc^Owe have 



(D k - V I-ifxI) l U k 



(k — r/) 2 + fi 2 



< 



nim(?; 2 ,(l + r ? ) 2 )(l + |//|) 2 
and it is not difficult to see that a similar estimate holds for \(Dq — nl — ijj,I)Uo\ 2 ■ Hence 



U k \ 



D — nl — i/il 



U 



l 2 xl 2 xl 2 xl 2 



< 



C 2 



A-e: 



4 Dichotomies near the core 



□ 



The system (7) can be abbreviated and written as 

V = (A-+B(s;\,p))U, 

where 



/ 



B(s;X,p) :=e 2 



(11) 



(12) 



0^ 

10 



\\-6(s)-p(s,-) 0/ 

We will show that the system (11) has an exponential dichotomy on the interval </_ = (—00, si]. 

To show this, we would like to apply Theorem 1 of [14]. This is not possible, however, since p £ 1Z is not smooth 
enough in s. We are interested in p small and consider therefore first p = 0, and show that the A-perturbcd 
system 

U' = (A_ +B(s;X,0))U (13) 
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possesses an exponential dichotomy in X on J_. Then we will use the implicit function theorem to show that 
also the system (11) possesses an exponential dichotomy. Note that from its definition, it follows immediately 
that B(s;X,0) e C(X). 

As 6 does not depend on <p, the system (13) decouples in Fourier space, just as the limiting system (8). Using 
the same notation in the Fourier spaces as before, we get for fceZ 



U' k (s)= A_(k) + B(s;\,0) U k {s). 



(14) 



As A_(k) = M k D k M k 1 , we rescale both U k and s to get estimates which are uniform in k. For k ^ 0, define 
r = \k\{s - si) and V k (r) = M^U^t/^ + Si). Then (14) becomes 



dr 



V k = 



D 1 + - n M k 1 B(T/\k\;X,0)M k 



V k 



(15) 



A short calculation shows that 



|Jl^ 1 fl(r/|fc|;A,0)M fc | < '—^ sup{l, |A - 9{s)\/k 2 }. 

* s< Sl 

Hence there exists a constant C such that for all k ^ , \M k 1 B(r/\k\; A, 0)M k \ < 2Ce 2T/|fe| and 

ro 1 r ° 2Ce 2r /l' £ l 



jT ^|Af- 1 S(r/|fc|;A,0)Af fc |dr<^ 



\k\ 



dr = C. 



(16) 



By the proof of the roughness theorem for ordinary dichotomies (see [6] for details), the system (15) has an 
exponential dichotomy which we denote by *5> k s (T,o), with constants K, n u = 1, k s = —1. We choose the 
dichotomy in such a way that Ran^l^i, s\; A, 0) C span{ei,e2}, where ej, j = 1,...,4 are the standard basis 
vectors of C 4 (again see [6]). This implies that the stable and unstable solutions satisfy 

\n(^r)V k \ < Ke-^- T ) \V k \, r < a < 0; 
\n(v,T)V k \ <Ke^-^\V k l a<r<0. 

The norm in X induces a norm on the Fourier space X k with 

\\U k \\ 2 Xk := ||^e ifc -|| 2 x = (k 2 + l) 2 ([^]i) 2 + (k 2 + l){[U k ] 2 ) 2 + (fc 2 + l)([U k } 3 ) 2 + {[U k ],) 2 . 

As seen in the proof of Lemma 3, M k is a linear homeomorphism between C and X k . Thus if we denote the 
exponential dichotomy of the unsealed system (14) by $^ s , then $^ s (s, t) = M k ^> k ^ s (\k\(s — si), \k\{t — si^M^ 1 , 
and they satisfy for U k £ X k 



\\§ s k {s,t)U k \\ Xk <Xe-W< s - f ' \\U k \\ Xk < K e-^\\U k \\ Xk , t < s < a i; 
||*2(M)£fclk < # e^"') ||C/ fe || Xfc < Xe^)||C/ fc || Xfc , s < t < 8l . 



(17) 



for some constant K 7 which is independent of k. 

For the central space, corresponding to k = 0, the scaling Vq(s) 



= ±es 



Uo(s) and the integrability of B(s; A, 0) 



shows that, for any e > 



|*o(M)fo| <Ke'^\U \, 
|*o(M)fo| <Ke-^\U \, 



t < s < si; 
s < i < si- 



(18) 



Thus for the full solutions, we can define the stable and center-unstable solutions 

$ s _(s,t;\,0)U = Y, m^t)Uke lk -, * < t < *i, 

fcez\{o} 

^ n(s,t)U k e ik \ t<s< Sl , 



$ cu (s,t;\,0)U = $o(s,t)U 
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and the unstable and center-stable solutions 

$ c _ s (s,t;\,0)U = $ o (s,t)U o + $fe(M)£4e ifc \ a < t < s u 

fcez\{o} 

A, 0)17= $t(s,t)U k e ik - , t < s < Sl . 

fcez\{o} 

These solutions are related to dichotomies for (13) in X on J_. 

Lemma 5. Lei — 1 = n s < k cu < < k cs < n u = 1 and A G R. T/ien </ie system (13) ftas an exponen- 
tial dichotomy in X on J_ wit/i constant K and rates n cu and n s , and another with constant K and rates k w 
and n cs . The dichotomies can be chosen such that Ran$l (si, si; A, 0) = P s and Ran $1 s (si, si; A, 0) = P cs . 
Moreover, for any t G (— oo,si] and Uq G X, the solutions $™(-, t; A, 0)[/o and $1 (•, i; A, 0)C/o belong to 
C°°((~oo,t);X) and C°°((t,Si);X), respectively. Similarly, the solutions <&"(■,*; A, 0)Uo and <f>ff(-, i; A, 0)U 
belong to C°°((— oo, t); AT) and C°°((f, Si); JQ, respectively. All solutions also depend smoothly on the parameter 
A. 

Proof. The scaling e ±r,s U for < 77 < 1 and the dichotomy estimates in (17) and (18) immediately prove the 
first part of the Lemma. The dichotomies satisfy Ran$l(si, s±; A),0) = P s and Ran$L s (si, si; A, 0) = P c + P s 
since we have chosen the Vf^ above to satisfy Ran\E'l(si, s±; A, 0) C span{ei,e2} (cf. the definition of D^). 

The smoothness with respect to s follows since 9 is smooth in s and smoothness in A can be proved using an 
implicit function theorem argument. First observe that for any A, A close to each other, the solutions $ c " and 
4> s satisfy the integral equations 



= - $™( s , t; A, 0) + $™(s, t- A, 0)) + (A - A) 
ft 

$™(s, r; A, 0)e 2r £ $ ra (r, t; A, 0) dr 
= - (a, f ; A, 0) + $1 (s, «; A, 0)) - (A - A) 
$1 (s, r; A, 0)e 2T J B $l (r, t; A, 0) dr 4 



$1 (s, r; A, 0)e 2r B $l tl (r, t; A, 0) dr 
$™(s, r; A, 0)e 2T B $l (r, i; A, 0) dr 
$1 (s, r; A, 0)e 2T S $ OT (r, i; A, 0) dr 
r; A, 0)e 2r £ <r (r, i; A, 0) dr 



for s < t < si and i < s < s\, respectively, where Bq is the matrix 

/0 OA 



Vl 0/ 



B n := 



(19) 



Define the function spaces 

X s := {<F 5 ; $ s (s,t) G £(X) is defined and continuous for t < s < si 
with ||$ s || s := sup e- KS ( s -*)||$ s (M)|| £(x) }; 

X™ := {$™ ; $ c "(s,f) G £(X) is defined and continuous for s < t < si 



with ||$ c 



c „ := sup e 

s<t<si 



™(- t )ll$™( s ,t)ll £(x) }. 



For A fixed, the integral equations can be written as F($ cu , $ s ; A) = 0, where F : X cu x X s x M -> X c " x X s . 
The estimates of the exponential dichotomies immediately give that F is indeed a mapping between those spaces, 
for example, 

^2 g 2si g K c "(s-t) 



$1 (s, r; A, 0)e 2T B $_ (r, i; A, 0) dr 



£(A') 



2 + k cu - K s 



The other integrals can be estimated in a similar way. Since Dr$cu t §s)F(§> cu (8, t; A,0), 3> s (s,i; A,0); A) = I, the 
implicit function theorem can be applied and the smoothness with respect to A follows immediately. □ 
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Remark 1. The ip -independence of 8 is not essential in Lemma 5. The lemma can be proved for smooth in- 
dependent functions 9 by using Theorem 1 of [14] and verifying the conditions (HI), (H2), (H3) and (H5) of 
that paper. 

Next, we prove four technical lemmas needed in the proof of the existence of exponential dichotomies for the full 
system (11) with p G 1Z. We work in exponentially weighted spaces, and for an unbounded interval Jc L and 
r\ G R , we let C V (J; X) be the space defined by 

C V (J;X) := {U G C(J;X); \\U\\c v ~ supe" s ||(7( S )|| A - < co}. 

s£j 

Hence Cq{J,X) is the space of continuous functions with an X-norm that is uniformly bounded in J. 
Lemma 6. Let J C J_ and pick u G C ( J; H 2 ^ 1 )) and p G L 2 (J; ff 1 / 2 ^ 1 )), then pu G L 2 (J; H 1 ' 2 ^ 1 )). 

Proof. We need to prove that for s fixed, 

\\p(s)u{s)\\ H1/ 2 {s i ) < C\\u{s)\\ H 2 (s i ) \\p(s)\\ H1/ 2 {S i ) . (20) 
Indeed, if this is proved, the claim follows, since 

11^11^(7^1/2(51)) = J \\p{s)u(s)\\ 2 Hl/2{sl) ds 

< C 2 ||p(s)||^i/2 (s i)||w||ff2( S i ) ds 

-f2 „„„ ||„,/'„M|2 / ll„/'„\l|2 



<C sup||u(.s)|| ff2(sl) JJ\p(s)\\ H i /2 ^) ds 

= C 2 |I u IIco(J;H 2 (S 1 ))II/ 3 IIl 2 (,/;H1/2(s1))- 

To prove (20), let u G i? 2 (5 1 ) and p G i? 1 ' 2 (S' 1 ) (we suppress the variable s for simplicity of notation). Let pt 
and Ufc be the Fourier coefficients of p and u, respectively. We have 

\\uf H2 = J2nl(l + k 2 ) 2 , 

iipii^/2 = E^( 1 + fc2 ) 1/2 - 



Then (up) k = £ , 6Z u jPk-ji an( l so II m pI| 2 ^i/2 — Sfcez (Ejgz u jPk-j^J (1 + fc 2 ) 1 ^ 2 . Let u and cr be the functions 
with Fourier coefficients Ufe(l + A: 2 ) 1 / 4 and pfc(l + fc 2 ) 1 / 4 , respectively. Note that v G _ff 3 / 2 (S' 1 ) and cr G L 2 (S 1 ). 

Now observe that 

1 + fc 2 = 1 + ((k - j) + j) 2 < 2(1 + j 2 ) + 2(1 + (k- j) 2 ), 

and hence 

(1 + fc 2 ) 1 / 4 < 2 1 / 4 ((1 + j 2 ) 1 / 4 + (1 + (jfc - j) 2 ) 1 / 4 ) 

for any j G Z. Thus 

iMiw < v^We^-^ 1 + j 2 ) 1/4 +Ew-i(i + ^ - j) 2 ) 1/4 ) 2 

fcez S'eZ jez ' 

< 2 ^x;((e + /) i/4 ) 2 + (e%^(! + ( fc - j) 2 ) i/4 X) 

= 2a/2(|MI| 2 + IIHl!*) 

<2>/2(sup |«M| 2 ||Hli=.+ sup |uM| 2 ||cr||l 2 ) 

^^(llfll^ibiiia + iiuii^ikiiio 

^c'dhii^ibiiw) 

for some constant C > 0. This completes the proof. □ 
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For each p G 1Z and s£ J_, let 

5B(s; ~p) := B(s; A, p) - B(s; A, 0) - -e 2s p(.s) B , 

where Bq has been defined in (19). Note that for any s G J_, 

II^B^PjlUtx) < sup ||e 2s p(s) m llia^i) < C sup sup |wi(y)|||e 2s p(s)|| L 2 (s i) 

II«iIIh2=i ll«ill„2=i (21) 

<C e 2s ||p(.s)|| i2(sl) , 

where we use the notation C for the different constants occurring. It follows that 



\\5B{s;p)f c{x) ds<C 2 r e 2s \\p( S )\\l 2{sl) e 2s ds<C 2 e 2 ^ / ||,5( S )|| 2 2(sl) e 2 * ds 

J — OO 

<cV^||pII|. 



(22) 



Lemma 7. For r\ G ( — l,« c "), where k cu is as in Lemma 5, pick UZ U € C V (J-; X) , and p G 7Z. Let s G J_. 
TTien i/ie integral 

I:= ( A-e A - p3( ~ s -^P s SB(T;p)Ul u (T)dt 

J —oo 

belongs to X . 

Proof. Let ff(r) := efo- x > T &B(<r; p)U± u (t). By the definition oi6B(r,p), 

H{t) = (0,0,0, -e(" +1 ^(r)u(r)) T =: (0,0,0,Mr)) r , 

where u(t) is the first component of UZ u (t). Then e r ' u G C ( J_; H 2 (S 1 )) and e'/5 G L 2 ( J_; i? 1 / 2 (5 1 )), and so 
by Lemma 6, ft. € L 2 (J_; H 1 / 2 (S 1 )). For fc G Z, let Hk(r) and /ifc(r) be the Fourier coefficients of -ff(r) and ft(r), 
respectively. Let P| := M k ~ 1 P^Mk- To show that / exists in X, it suffices to show that {Ik}kei, & X (see (9)), 
where 

/ fc := f e( 1 -") r M Jfc D fc e D *^f a -' r )M i 7 1 P'fr fc (T)dT 



^ f e( 1 -") T e-l fc K- T ^ fc (T)dT(0,0,l,-|A;|) 



7' 



We therefore need to prove that 



— oo 



o^-V)r e -\k\(.s-r) hk ^ dT l £[ 2 

) fcez 

Using that rj < and ft. fe G L 2 (J-) (as h G L 2 (J_; i7 1/2 (S' 1 ))), we note that 

/s / /-s \ 1/2 / „ s x 1/2 

eCi-^e-lfcl^-r)^^)^ < [jfeje-lfcl- (J e 2(i-„+|H)r dT j Ij hk (jf dr ^ 

\/2(l -?7+ |fc|) 

^e^-^ci + lfcl 2 ) 174 !!^!!^.). 

Since {(1 + |fc| 2 ) 1//4 ||ftfc||L 2 (j_)}fcez G I 2 , the proof is complete. □ 

Lemma 8. for — 1 < r\ < k cu , where k cu < is as in Lemma 5, pick f/f" G C n (J~;X) and p <G 1Z. Then the 
integrals 

f $ s _(s 7 T;\,0)SB(T;p)U m (T)dT and f B(s; A,0)$!_(s,r; X,0)SB(t; p)U1 u (t) dr 

J —oo J —oo 

exist in X for each s G J- . 
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Proof. We use (21) and compute 

f ^ s _{s 1 T-\,0)SB(T;p)U r i u {T)dT 

J — OO 



X 



<C / ||$i( S) r;A,0)|U w e^||/5(r)|| L2(sl) ||C/™(r)|Udr 



< C I e 2T e KS(s - T) Wn 



\\p(r)\\msn\\Ut u (r)\\xdr 



<C\\UZ u \\c n (j^x)e K 
<C\\Ul u \\ Cv{J _. x y 



1/2 / r s \ 1/2 

e2( W„ J))rrfr) ( / e 2.^ (r)|| 2 2(si)dr 



— OO 

1 



" ll/5| 



R ■ 



y/2{l - K° - rj) 

Using that B(s; A, 0) <G C(X), it follows that both integrals converge in X. 



□ 



Lemma 9. Let —1 < r/ < k cu , where k cu < is as in Lemma 5. Let UZ U G C n (J-; X) and p £ 1Z. For every 
s G J- , the integral 

f A-$l{s,T;\,0)SB{T;p)U cm (T)dT (23) 

J — OO 

exists in X . 

Proof. By (3.1) of [14], for r < s < Sx, 

$l(s,r;A,0) = e A - ps(s - T) P s - f e A ~ F "( s ~^ P S B(£; A, 0)$™(^, r; A, 0) d£, 

J — OO 

+ £ e^^-^P^; A, 0)*i(& r; A, 0) d£ 
- f 1 e A - p " u ( s -®P cu B(£; A, 0)$i (f , r; A, 0) df. 

J s 

Note that we used here that Ran$l (si, si; A, 0) has been chosen so that it coincides with RanP s . 

Substituting this into (23), we have four integrals to estimate, the first of which was dealt with in Lemma 7. 
The other three integrals are 

h~ f A- T e A - ps ^P s B(^X,0W^(^T;\,0)d^B(T;p)U c _ u (T)dr, 

J — OO J — OO 

h~ f A_ f\ A - pB ^P s B{^\ 1 ())^,T-\ 1 ())d^B{r-rp)U^{T)dT, 

J — OO J T 

h-= f A_ r e A - p "'^P cu B(^X 7 0)^t(^T;X,0)d^B(T-,p)Ul u (r)dT. 

J — OO J s 

We carry out the calculations for ii, since the others are similar. Let (<j>ji(£, r)), j, I = 1, . . . , 4, be the entries of 
the matrix corresponding to ( I>'i u (£, r; A, 0), and as in the proof of Lemma 7, let h{r) = —e^ r,+1 ' T p(r) uir). Recall 
that h G L 2 (J_; i/ 1/2 (5 1 )). A short calculation shows that 







B(£; A, 0)$l"(e, r; A, 0)SB(t; p)Ui u (r)e 



cur \ „2f+(l-r;)r 







Note that (^34 (^,r) and 014 (£,r) map L 2 (S 1 ) boundedly into H 1 (S 1 ) and 7J 2 (5' 1 ), respectively, and that by 
Lemma 5 for £ < r < si 

\\MZ,t)\\c (l *.hi) < Ke" *^, 
\\M£,t)\\c(l*;h>) <Ke^^ T K 
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Introducing the notation /(£,r) := e- K "'^ T ^(f> 34 (^ t)1i(t), and ff(£,r) : e^"^) (A - 0(r))0i 4 (£, r)ft(r), we 
note that max(||/(£,r)||#i, ||<7(£,t)||#2) < i : f||/i(T)|| i 2, for £ < r < si. The Fourier coefficients of /(£, r) and 
<?(£,r) are denoted by t) and <7fc(£, r), respectively. 

To prove that Ji S X, it suffices to prove that { Jfcjfcgz € X, where 



Jk 



M k D k 



e 25+(l-r,)r e K ct '(C-r) e I5 fc P=( ;i -5) A/ -lp S 



\ 



V&(£,t)/ 



e 2C+(l-r,)r e K<="(C-r) e -|fe|(s-0 



oo </ — oo 



?fc(C,T) 

V-|fc|?*«,r)y 



d£ dr, 



where P£ = M k 1 P k s M k as before. The first component of J k can be written 
1 



-|fc|« 



= (1_ k <="_ i))t 



3 (i+« c "+|fc|)£ p ? 



" 2 



1 



< 



1 



P (2-»?)s 



1*1) 



1/2 



4 ^(1 + ^" + |fc|)(2+ |fc| -7?) 



G>0 •/ — OO 



1/2 



The square of the l\ norm of the first component of {Jk\k& can then be estimated by 



2(2- V )s 



16 H, ( 

fce 7 v 



^ (1 + K c 



1 + fc 2 



OO </ — oo 



feez 



16(1 + K cu ) 
e 2(2- V )s 

32(1 + k cv 



h|fc|)(2+|fc|-7 ? ) J_ 

e^h k { T fd^dT 



(l + k 2 )e 2 ^M,r) 2 d^dT 



oo J — oo 



p 2(3- 9 j)a I 
<^- rl|^ 2 



^ /' e 2 ^ fc (r) 2 dr 



32(1 + k c1 



(24) 



The other three components of {Jk\k& are estimated in a completely similar way, and by adding these estimates 
we see that T\ £ X. □ 

We are now ready to prove the existence of exponential dichotomies for the full system. 

Theorem 2. Let —1 < k s < k cu < and < k cs < k u < 1. Then there exists a neighbourhood U 
of in 1Z such that for any p £ Li and any A £ M, the system (5) /ias an exponential dichotomy on J_ 
TOi/l constants K and rates k cu , k s , and another with constants K and rates K u , n cs . Moreover, the pro- 
jections and evolution operators depend smoothly on A £ R and p £ hi. The dichotomies are denoted by 
$f_(s, t; A, p), $™(s, t; A, p), and $l s (s, t; A, p), (s, t; A, p), respectively. The associated projections will be 
denoted by P0_(s; A, p)(:= & s _(s,s; A, p)), P^ u (s;X,p), Pl s (s; A, p), and P"(s; A, p) , respectively. 



Proof. We will show that there exists a neighbourhood of in 1Z such that if p belongs to this neighbourhood 
then there exist exponential dichotomies for the system (11) with this p. Let Uq £ X and t £ J_ be fixed but 
arbitrary. We will use the implicit function theorem to solve the system of integral equations for the pair of 
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(25) 



functions (£7™, Ut) as functions of the parameters Ael and p € 1Z near 

= <S> au (s,t;\,0)U -U° u (s) + f $!_(s,r; X,0)SB(t; /5)[/ ot (t) dr 

J — OO 

-J ^ m { Sl T-\Q)5B(T-p)U m {T)dT 
+ J <£>™( Sl T;\,0)6B(T;p)Ul(T)dT, for s < t < s 1; 
= &L(8,t;\,0)U -UL(a)- f $l(s,r; A,0)<55(r; p)U1 u (t) dr 

J — OO 

+ ^ $l(.s,r;A,0)<5B(r;p)[/l(r)dr 

-f $™(s,T;\,0)5B(T;p,0)UL(T)dT, for t < s < Sl . 

J s 

By Lemma 5, the dichotomies <J>f"(s, £; A, 0)C/o and $i(s, A, 0)£/o exist and have constants K, k s = —1 and 

k cu e (-1,0). 

Let 77 € (k s ,k cu ) and rewrite equation (25) as F(Ui u ,Ui; A, p) = 0, where F : C v ((-oo,t}; X) x C v ([t, Si]; X) x 
ixK-> C,((-oo,i];X) x C v ([t, si]; X) is the right hand side of (25). 

We first verify that F is indeed a map between the above spaces. We do the estimates for the first integral 
in the first equation of (25). The other estimates are similar. Lemma 5 gives that for any s £ (— 00, t] and 
U cu G C v ((-oo,t];X): 

f <S> s _(s,T;\,0)5B(T;p)Ul u (T)dT 

■/-oo X 

<K sup (e^\\U c _ u (r)\\ x ) f e^-^ s -^\\dB(T;p)\\c(X)dT 

3,3] J —OO 

CK-oct],*) f ( e 2(KS -" )(s - T) + ll^(r;/5)||^ w ) dr 



t£ ( — 00, sl 



(26) 



< ^l|c/ ra !lc, j( (-oo, t ],x) 



2(7? -K*) 



e 2si H/5||| 



where we have used (22). After taking the suprcmum over all s G (—00, si] we sec that the function defined by 
the first integral in (25) belongs to C V (J-, X). Using similar estimates for the other integrals, we can conclude 
that F is indeed a map between the spaces as stated. 

That F is smooth with respect to A and p follows since the evolution operators ( &'i tl (-, t; A, 0)Uq and $!_(•, t; A, 0)Uo 
are smooth in A by Lemma 5 (using that the H 1 norm is weaker than the C 1 norm on bounded intervals), and 
since 5B depends smoothly on p (indeed, 5B is a bounded linear mapping with respect to p). Note that 

F(*«(-, t; A, 0)U , $1 (•, t; A, 0)U ; A, 0) = 0. 

The Frcchet derivative of F with respect to its two first variables evaluated at ($l tl (-, t; A, 0)Uq, i; A, 0)[/o; A, 0) 
is —I on Crj((— 00, i]; X) x C^([t, si]; X). In particular, this derivative is a linear homeomorphism on this 
space, and so the implicit function theorem is applicable, and we obtain solutions $!"(-, t; A, p)Uo := f/f" and 
$L( - , t; A, p)Uq := Ut of the integral equation (25), which exist in a neighbourhood of (Ao, 0) in R x TZ. Smooth- 
ness of these solutions with respect to parameters also follows from a corollary of the implicit function theorem 
(see e.g. [4, p.115]). 

Next, we need to verify that $™(-, t; Ao, p)Uo and $f_(-,t; Xo,p)Uo are weak solutions of (11), and that they 
satisfy the conditions of Definition 2. Wc first check that $™(-, t; A, p)Uo is a weak solution on the interval 
(— 00, i]. By Lemma 5, $™(-, t; A, 0)Uo is a C°° solution of 

U' = (A- + B(s;X,0))U 
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on (— oo,i], and hence it is also a weak solution of this equation. Next we deal with the integral terms. For the 
first integral we use the abbreviation 

g(s):= f f(s,T)dr, with f(s, r) = *! (a, r; A, 0)6B(t; p)*~(t, t; A , p)U . 

J — oo 

Thus / is C°° in the first variable and L 1 in the second. From its definition, it follows immediately that g is 
continuous. We will see that g is weakly differentiable and that 

g'(s) = f(s,s)+ f ?f(s,T)dr. (27) 

J-oo OS 

In order to prove this, we need to check that the integral on the right hand side of (27) exists, and that the 
equality (27) holds. The integral in the right hand side of (27) is 

/ ^(s,r)dr = f (A. + B( S ;X,0))^ s _{s,T;X,0)6B(T;p)^(T,t;X,p)U o dT, 

J — oo OS J —oo 

and it exists in X by Lemma 8 and Lemma 9. 

Next, we calculate the distributional derivative of g and let V € oo, t];X) be a test function. Then by 

Fubini's Theorem and integration by parts 

g'(s)V{s)ds = - I g(s)V'(s)ds 

f(s,r) dTV'(s)ds 



oo J — oo 



oo 

t 



f{s,T)V'{s)dsdT 



f(r,r)V(T)+ I ^( s ,T)V(s)ds) dr 



d_ 

ds 



/l( /(s ' s) + f 00 % { ' %T)dT ) v{s)ds ' 

and we see that the weak derivative of g is indeed given by (27). Hence 
■ f <£ s _(s,T;\,Q)6B(T;p)4>™(T,t;\,p)U dT 

J — OO 

= (/ - Pl u (s; A, 0))SB(s; p)$1 u (s, t; A, p)U Q (28) 
+ / (v4_ +B(s;A,0))$l(s,r;A,0)(5B(r;p)$i u (r,t;A,p)C/odT. 

J — oo 

We have already noticed that g(s) = f* $1 (s, r; A, Q)SB(t; p)$™ (r, t; A, p)Uo dr is continuous, and so it belongs 
to Lf oc ((—oo,t];X). The right hand side of (28) also belongs to Lf oc ((—oo,t];X) since the first term belongs to 
if oc ((— oo,t];X) and the second term is continuous on (— oo, t\. This shows that 

f $1 (s, t; A, 0)6B(t; /5)$™(t, t; A, p)U dr 

J —oo 

belongs to Hl oc ((—oo,t];X). 

Similar calculations for the other integral terms of the first equation of (25) show that these are also weakly 
differentiable on (— oo,t] and belong to Hj- oc ((— oo, t], X). After adding the terms up, we conclude that 

^ c _ u (s, f; A, p)U = (A. + B(s; A, p))$™(s, t; A, p)U , 
ds 

i.e. $™(-,t; \,p)Uo is a weak solution of (11). 

Similar calculations for the terms of the second equation of (25) show that $!_(-, t; A, p)Uo is a weak solution of 
(11) on the interval [t, s±]. 
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Finally we check that the conditions of Definition 2 are satisfied. A similar computation as in (26) also shows that 
the estimates in (i) and (ii) of Definition 2 are satisfied for $™(s, t; A, p) and $1 (s, t; A, p) for any k cu and k s such 
that — 1 = k s < k s < k cu < k cu < 0. Since k cu can be taken arbitrarily close to 1, the same is true also for k cu . 
Note that (iii) of Definition 2 is satisfied with Pf"(s; \,p) := $ cu (s, s; A, p) and Pi(s;X,p) := $!_(s, s; A, p). □ 

To finish this section, we derive some more details about the solutions of (5) in the case when A = Ao and p = 0. 
We are particularly interested in the solutions on J_, and we study the exact growth/decay rate of solutions as 
s — > — co. As we have seen before, the space X decouples into a direct sum of 4-dimensional pairwise orthogonal 
Fourier subspaces Xk, and that since 9 is radially symmetric, the subspaces Xk are invariant both under the 
flow of (5) with p = and under the flow of the asymptotic system (8). 

Lemma 10. Let ej, j = 1, . . . ,4 be the standard basis of C 4 and consider the four- dimensional invariant central 
space corresponding to k = of the unperturbed equation obtained when p = and A = Ao in (5). Then there 
exist two unique solutions Uoj(s) with j = 1, 3 such that 

lim U ,j(s) = ej, j = 1,3. 

s— > — oo 

We may also pick two solutions Uqj with j = 2,4 which grow algebraically as s — > — oo and satisfy 

lim -U ,j{s) = e^, j = 2,4. 

s->-oo s 

The solutions Uqj, j = 1, ■ • ■ ,4, are linearly independent. 

Proof. It is straightforward to check the assertions of the lemma, using [5, Chapter 3.8]. □ 
In Section 6 we will specify the solutions {To, 2 and C/0,4 using the adjoint system. 

Lemma 11. For every k £ Z \ {0}, there exist solutions Ukj of (5) with p = and A = Ao such that (together 
with the solutions specified in Lemma 10 for k = 0) we have 

span{[/ fej (si); k G Z, j = 1, . . . , 4} = X, 

and for s — > —00, 

e W{s - Sl} U kA (s) -> (-l/fc 2 ,l/|fc|,0,0) T e*, 

e |fc|C— x) C^ >2 (49) — ► (0,0, — l) T e**", 
e-l fc K s -^)[/ M Gs)^(l/fc 2 ,l/|fc|,0,0) T e*, 
e -l*l(«-»0[/- M (a) _> (0,0,l/|fc|,l) T e lfc '. 

Proof. As seen in the beginning of this section, the system (5) with p = and A = Ao leaves the subspaces X k 
invariant. The estimates on the matrix B(s; Ao, 0) in (16) now show that there are solutions of (14) (and hence 
of (5)) which converge to the solutions of the system at infinity, see e.g. [5, Chapter 3.8]. In Section 3 we have 
seen that (8) has two solutions in X k with decay rate e^ s and two with growth rate e~' fc ' s for s — s- —00. A 
comparison with the eigenvectors of -A_(fc) in Section 3, we obtain solutions Uk,j, with k G Z\{0} andj = 1,...,4 
with the desired properties. □ 

Next, we perturb the solutions Uqj(s) with j = 1, 3 described in Lemma 10 to solutions of (5) for all sufficiently 
small potentials p. First, we will show that the four-dimensional central subspace corresponding to k = persists 
in (5) as the intersection of the ranges of $i s (si, s\; A, p) and $™(si, S\; A, p). Note that the difference between 
the operators A(s; A, p) and A(s; Ao, 0) in (5) is 

A(s; A, p) - A(s; A , 0) = r'(s) 2 (A - A - p)B Q e 2s (A — Ao - p)B 0) 

as r'(s) 2 = e 2s for s < si (see (19) for the definition of Bq). 
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The function e T p(r) belongs to L 2 ( J_, H 1 / 2 ). By Lemma 5, e T p(t)u\{t) also belongs to this space. Thus 
\\e T p(t)B U\\x € L 2 (J_) and \\e 2T p(t)B U\\x is the product of an L 2 function and the exponentially decaying 
function e T . 

This allows us to use the Gap Lemma as in [16, §4.3 and (4.12)] and [3, Proof of Lemma 4.1] to show that (5) 
has two linearly independent solutions U§ b j(s; A, p) for j = 1,3 that converge to e 3 as s — > — oo, and two other 
solutions which grow algebraically. In fact, the results in these works show that any linear combination of the 
bounded solutions U§j(s; A, p) with j = 1, 3 can be found as a fixed point of the equation 

U(s) = $ m ( S ;A ,0)(7 cb + f <S> c _ s ( S ,T;\ 07 0)e 2T (\-\ -~p(T))B U(T)dT (29) 




r; A , 0) e 2r (A - A - p(t))B U(t) dr 7 



where U§ b belongs to the unperturbed bounded central subspace spanned by Uqj(si) for j = 1,3. We denote 
the fixed point by U_ (s; A, p, U^ b ) and write 

Pi b ( Sl ;\pM b := Ul b (s 1 ;\,p,Ug>) = U£ b + H $ c J( Sl , r; A , 0) e 2r (A - A - p{t))B U?(t; A, p, U^ b ) dr. (30) 

J — oo 

Similarly, we can use (25) to describe the solutions of (5) with exponential decay as s — > — oo by 

P?( Sl ;\,p) = Pl I (.si;A o ,0)+ f ' $L a ( Sl) r; A , 0) e 2r (A - A - p(t))B *_ (t, s i; A, p) dr. (31) 

J — oo 

These results will be used later to characterize cigenfunctions of the perturbed operator. 



5 Dichotomies for the far field 

The method of Section 4 is not available for determining the dichotomies for s large. Going back to (2), we 
observe that 6 and p have support in a ball with radius r\, and thus for r > T\ the eigenvalue problem (2) reduces 
to (A 2 - X)u = 0, which can be factorized: 

(A - VA)(A + \fx)u = (A + VA)(A - \fx)u = 0. 

Expanding u(r, ip) as a Fourier series in the angular variable tp, we see that the Fourier coefficients u& satisfy the 
differential equations 

9 2 19 k 2 rr\ ( d 2 Id k 2 /-A _ „ 

TT + -^ 2 ~ v*) \ ITS + -* T + VA Sfe = 0. 32 

or z r or r z J \ or z r or r l J 

For k fixed, this is a fourth order linear ODE, so it has a four-dimensional space of solutions. The general solution 
can then be obtained as a linear combination of the solutions of 

\ dr 2 ^~ r dr r 2 k ^ ^ 

and the solutions of 

or z r or r z J 
so that the general solution of (32) is given by 

Sfc(r) = dhiX^r) + C 2 K k (X^ 4 r) + C 3 J k (X^ 4 r) + C 4 F fc (A 1/4 r), 

where J k and Y k are Bessel functions of the first and second kind, respectively, which satisfy equation (34), and 
Ik and Kk arc modified Bcsscl functions of the first and second kind, respectively, which satisfy equation (33). 
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For r > r3 = S3, we have s = r. Thus, we can define the systems corresponding to equations (33) and (34) with 
the variable s for s > S3 as 



1 (35) 
—~ + v A ui u 2 , 



fc 2 /-\ 1 (36) 
— - VA iti it 2 , 



and 



respectively. We will consider these systems for s > r% and define <j>k{s,t) and tpk(s,t) to be the evolution 
operators corresponding to the systems (35) and (36), respectively. After deriving dichotomies for those systems, 
we will derive dichotomies for the original system (5). 

To derive dichotomies for (35) and (36), we introduce for s > n the function spaces X s — H 1 (S 1 ) x L 2 (5' 1 ), 
Y s = H 2 ^ 1 ) x H^S 1 ) and Z s = H 3 ^ 1 ) x H 2 ^ 1 ) with norms 



s- 

- : = ^IKIIjj'cs 1 ) + IKIIiri(si) + ll u 2!lffi (S i), 



Nl|« := ^IMIi? 3 ^ 1 ) + IKIIiP(Si) + Il u 2!|ff2 (s i). 

We decompose the spaces X s , Y s , and Z s into their Fourier subspaces 



X s = ® keZ X s k , Y s = ® k ^Y k s , and Z s = ® keZ Z s k , (37) 

where 

X s k = Y k s = Zl = {(oe*, &e*) T 5 «, & € C}, 

and the completion in (37) is in the respective norms of X s , Y s and Z s . The norms on X|, F fc s and Z| are given 
by the restriction of the norms of X s , Y s , and Z s respectively, and so 

||(ae* te«ri&.= (i + 5) H 2 + H 2 ; 

||(ae*,6e ife -) T ||| a (l + fc 2 ) ||(ae ife \ 6e*) T ||| 8 ; ( 38 ) 

k k 

||(ae*,6e*) r ||| s (l + fc 2 ) 2 H(^ fc -.^ fe -) T ||| s . 



For each e G (0,A 4 / 4 ), we now prove the existence of a time-dependent exponential dichotomy for (35) with 
constant K > and rates k s = —(A 1 / 4 — e) and k u = (A 1 / 4 — e) that are independent of fc. For (36), we will 
show that the evolution operator always acts in the center-unstable manifold and derive that its growth can be 
bounded by any exponential. 

Lemma 12. There exists an Cq > such that for any e G (0, eo) there exists a K > such that for any k G Z 
and A G (Ao/2,2Ao) i/iere exists a time- dependent exponential dichotomy of (35) on J + so £/ia£ 0|(s,t;A) and 
(j> k {s, t; A) satisfy 



H s k (s,t;X)\\L(xt,x t) = 1101(^)11^,^) < xe-^ 174 -^-') s > t > n, 

n u k (s,t;\)\\ c( x ljtt) = HUs,t;X)\\ ci Y*,Y°) < ^- (Al/4 - £)(! " s) t > s> r x . 



Proof. Let (ui,U2) T satisfy equation (35). To get estimates which are uniform in fc, we follow [18], and let 

U\ 1/2 



/ fc 2 \ 

h(s) := ( VA+ — j ui(s). 
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Note that 



m(l,v^V2)||(«i,U2) T ||jt. < || u 2 ) T \\c- < max(l,y2A^)||( Ul , U2 ) T ||^ 



and that the constants above are independent of k and A. This shows that, when using the new variables u\ and 
U2, we can use the standard norm in C 2 . 



Next, we rewrite the system (35) in the new variables u±, u 2 '. 



21 = V\ 



2\ 1/2 



U2 



1 

-~U 2 



rr k 

v A + — 



•S" 

2\ 1/2 



fc 2 ( rr k 2 
v A H — - 

s z 



Ui, 



Ux, 



Now, we change the independent variable by making the substitution dr/ds = (VA + fc 2 /s 2 ) 1 / 2 . We write s(r) 
to describe the dependence of s on r. We then obtain (where ' now denotes differentiation with respect to t) 



k 2 



s(r) a 



1 



Ul 



s(r) 



Va + 
Va 



k 2 



S (r) 2 
k 2 



-3/2 



*(t) 5 



-1/2 

u 2 . 



(39) 



Noting that s(r) — > oo as r — > oo we find that the limiting system at +oo is 

u[ = u 2 , 

Uo = Ux, 



(40) 



which is independent of fc. The matrix associated with this system has eigenvalues ±1. Hence equation (40) 
possesses exponential dichotomies with n u = — n s = 1. To get estimates for the perturbated system (39), we will 
use the estimates 



s(r) 3 
1 



ST 



a/A 

Va 



s(r) 2 

fc 2 
s(r) 2 



-3/2 



-1/2 



J. 1 fc 2 /s(r) 2 1_ 



< 



This estimate is uniform in A in a neighbourhood of Ao- The roughness theorem for exponential dichotomies [6, 
Chapter 4] now guarantees the existence of an exponential dichotomy also for the system (39), and we denote 
the corresponding evolution operators by (f>l(cr, r) and <^(o\ t). For each positive e sufficiently small, there exists 
a K > such that 

||^M|U (C 2) < Jf e -(l-e)(<r-r) ff > Tj 

||^(^t)|| £(C2) < # e -(i-?)(-') a < T . 
Moreover, X does not depend on A in a neighbourhood of Ao or on k E Z. 

It remains to translate this result back to the s variable. We write s = s(a) and t = s(r). Note that ds/dr < A 1 / 4 , 
and so by the chain rule we have for s > t 

UlM\\c { Xi,X t ) < CU s k (a,r)\\ c(c2) < /fe-Ci-^C-r) 

< tfe-a-^C-t) < ife"( Al/ ^ £ )( s -*), 

where we have put e = eA 1 / 4 . A similar calculation proves that for t > s 



\4>t(s,t)\\c { xi,x») <Ke 



-(A 1//4 — e)(t — s) 



The estimates for Yj° also follow from these estimates, since it is only a matter of multiplying both sides of the 
inequalities by a factor (1 + fc 2 ). □ 
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Lemma 13. Let e > be given. Then there exists a K > such that for any k G Z and A G (Ao/2, 2Ao) we /iawe 
||V>fcM;A)|| £( ^ ii:|) =||^fc(*,t;A)|| jC(52i? . ) KKe 6 ^ for t>s> n . 

Proof. First we note a scaling invariance in (36). If (ui(s, i), m 2 (s, i)) is a solution of (36) with A = 1, then 
(ui(s,t),U2{s,t)) = (ui(A^ 4 s, \\^t), A^ 4 U2(A^ 4 s, \\^t)) is a solution of (36) with A = Ai. So it is sufficient to 
prove the estimate in the Lemma in case A = 1. Using the explicit expressions for the solution in terms of Bcsscl 
function, it follows that, for A = 1, ip k (s,t) is given by 

i>k(s,t) 



(J k {s) 


Y k (s)\ 


J k (t) 


Y k (t)\ 






vm 





= Trt J k (s) Y k (s)\ ( Yl(t) -Y k (t)\ 
' 2 Yl{s)j \-J' k {t) J k {t) j 

_ Trt f a k (s,t) b k (s,t)\ 
~ 2 \c k {s,t) d k {s,t))' 

where 

a k (s,t) := (J fc ( fl )y^(*)-y fc ( S )4(t)), 
b k {s,t) := -J k {s)Y k {t)+Y k (s)J k {t), 
c k {s,t) := (J' k ( S )Yl(t)-Yl(s)J' k (t)), 
d k (s,t) := (-J' k (s)Y k (t)+Yl(s)J k (t)), 

and we have used that the Wronskian of J k (t) and Y k {t) is ^ [1, (9.1.16)]. Writing u = (ui,it 2 ) T G X kl we have 

||V>fc(Mkll~s - — ( U + J (afc(s,t)Mi + 6fc(s,i)w2) 2 + (cfc(s,t)wi+d fc (s,i)w 2 ) 2 
and so we need to show that there exists a K > such that for every u\ , «2 G K, A; G Z and t > s > n , 
+ MMW+MM)w2) 2 + (c fe (MK+4(M)«2) 2 ) <K 2 e 2e{t ~ s) (Jl + ^\ u\ + i 

By choosing u\ and ui appropriately, we note that this inequality holds if and only if the following two inequalities 
hold for some K > and all k G Z and i > s > r±: 

+ J) « fe (M) 2 +c fc (M) 2 ) i 2 < tfV*"') + 

((l + 5) MM) 2 + ^(M) 2 ) * 2 < ^ 2 e 2e(t " s) . 

To simplify further, we note that the above two inequalities hold if there exists a constant K > such that for 
t>s>r u 

U2 I ~V2 

l + -\a k (s,t)\t<Ke^ s ^l + ¥ , 



k 2 

l + -\b k (s 7 t)\t<Ke^- s \ 

\c k (s,t)\t<Ke< t - s ^l+ 1 ^, 
\d k (s 7 t)\t<Ke^- s \ 

First we will prove the fourth inequality of (41). Let 

f k (s,t) :=e-^~ s hd k (s,t). 

We need to show that f k (s,t) is uniformly bounded for k G Z and t > s > r±. Since J- k = (—l) k J k and 
Y~_fc = (— l) k Y k , it is sufficient to consider k G N. 
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We start with showing that f k is bounded for k <E N fixed. In the slightly smaller sector t > (1 + S)s > (1 + S)ri 
(where 5 > is arbitrary), we have \f k (s,t)\ — > as s 2 + t 2 — > oo, or equivalently, as t — > oo. Indeed, Y k (s), 
Jfc(s), V^(s) and J' k {s) are bounded by a constant Cfc for s > n [1, (9.2.1)], and so 

|/fc(M)| < Gle-< l - s h < C 2 k te-^ t/(1+s) -> (42) 

as i — > oo. Furthermore, y / sifc(s), >/sJk(s), y/sY k (s) and y/sJ' k (s) are bounded by a constant for s > ri [1, 
(9.2.1)], and so for n < s < t < (1 + S)s we have 

\fk(s,t)\ < e-^- s ^(l + S)std k (s : t) < VTTSD 2 . 

Altogether this implies that fk(s,t) is bounded in the whole sector t > s > n by a constant, possibly depending 
on fc. 

To show that in fact f k {s,t) is bounded by a fc- independent constant, we consider s > T\ as being fixed for the 
moment. First we note that in (42), we proved that for fixed s and k, fk(s,t) — > as t — > oo. Thus for s > r± 
fixed, the function f k (s,t) attains its maximum in an interior point t > s or at the boundary t = s. We use a 
method by L. Landau [11] to analyze the behaviour of fk(s, ■) at its critical points. At a critical point, we have 
dfk/dt(s,t) = 0. By equation (11) of [11], at the points where df k /dt = we have 

hA(.,o*-» ; ^ ? |(.-^)^w). 

where Ak(t) = J °° K (2t coshr)e~ 2kr dr, and Kq is a modified Besscl function of the second kind, satisfying 

K (x) = / e~ xcoshT dT. 
Jo 

In particular (since 2tf k (s,t) 2 /e 2 ^ t ~ s H 2 > 0), f k (s,t) 2 is decreasing in /c at a point where dfk/dt = if and only 
if e~ 2e ( t-s )£ 2 j4fc(t) is decreasing in t. Note that A k (t) is monotonically decreasing for t > 0, and since e - 2e (*- s )£ 2 
is monotonically decreasing for t > 1/e, we conclude that |/fc(s, -)| is monotonically decreasing in fc at its critical 
points for t > 1/e. 

Note that in the case when r\ > 1/e, we have proved that if the maximum of fk(s, •) occurs for t > s > ri, then 
the maximum is decreasing in fc, and hence stays bounded as k increases. At the boundary t = s > ri, we have 
fk(s,s) = 2/-7T, which is independent of fc. As each function f k (s,t) is bounded, in particular /i(s,t) is bounded, 
it follows that fk(s, t) is bounded in the whole sector t > s > r\ by a /c-independent constant for all fc £ Z. This 
shows that is uniformly bounded in the case when r\ > 1/e. 

When r\ < 1/e, we also need to estimate fk(s,t) in the triangle 1/e > t > s > r\. Here we use the estimate 
e -e{t-s) t < ij e It f n ows t h at \f k ( s ,t)\ < \g k (s,t)\, where 

. 9fc (.s,i) := -(Yl(s)J k (t) - Ji(a)n(t)). 
e 

Applying Section 3 of [11] wc conclude that g k (s,t) 2 is decreasing in k at the points where dg k (s,t)/dt = 0. 
Furthermore, gk(s, t) — > 0, for t — > oo and gk(s, s) = < for s > ri. The proof of the fourth inequality of 
(41) is complete. 

Next, we prove the second inequality of (41). By [1, (9.1.27)] we have 

(l+ *) Ms) = Y k (s) + i(y fc _i(s) + F fe+1 ( S )), 
(l + ^ Jk(s) = Ms) + \{Jk-i{s) + J fe+ i(s)). 

Note that 

( k 2 \ 1/2 k k{ k 2 \ 1/2 
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and so the second inequality of (41) is equivalent to 



< Ke e{t - S \ 



t ({Y k {s) + il*_x( S ) + iy*+i00) Mt) ~ (Jk(a) + \jk-i{s) + ±J k+1 (s))Y k (t)\ 

To prove the inequality, we use the same method as above, with 

f k (s,t) := te-^-^ (^Y k (s) + + ±Y k+1 (s))J k (t) - (J k (s) + ij fc _i(«) + \ J k+1 (s))Y k (t) 

and 

g k (s,t) := tz- s& — ({Y k (s) + ^Y^s) + ±Y k+1 (s)) J k (t) - (j k (s) + ±J k -x(s) + \ J k+ i(s))Y k (t)) . 

Then \ f k {s, t)\ — > as s 2 + t 2 — > oo in the sector t > s > r\ and as above, f k {s, t) 2 is decreasing in k at the points 
where df k /dt = if t > 1/e. In the triangle 1/e > t > s > r\, \f k {s,t)\ < \g k (s,t)\ and g k is decreasing in k at 
the points where dg k /dt = 0. Finally, at the boundary points where s = t > r\, we have f k (s, s) = g k (s, s) = 0. 
We conclude that the second inequality is valid in the whole sector t > s > n. 

For the third inequality of (41), we use the last identity of [1, (9.1.27)] which shows that the inequality is 
equivalent to the two inequalities 

k ' ^ 1/2 

t 



''■\J'k(s)Y k (t) - Yl(s)J k (t)\t < Ke<^ (l + ^) , 
\J' k {s)Y k+l {t) Y k \s)J k+l {t)\t < Ke^ (l+ ^ 



The first of these inequalities follows from the fourth inequality of (41), and the second can be handled as in the 
proof of the second and fourth inequalities after noting that at the boundary where t = s > n we have 

k 9k 
\J' k (s)Y k+1 (s) - Yl(s)J k+1 (s)\ = -\J' k (s)Y k (s) - Yl(s)J k (s)\ = — , 

thus f k (s,s) = — r=^ < -• We omit the details. 

It remains to prove the first inequality of (41). This can be handled as the second inequality by using 

J k (s)Y k \t) ~ Y k {s)J' k {t) = *(J fc (s)n(t) - Y k (s)J k (t)) + (Y k (s)J k+1 (t) - J k (s)Y k+1 (t)) 

and splitting the inequality up into the two inequalities 

/ k 2 \ 1/2 k ( k 2 \ 1/2 

(1 + - t \J k { S )Y k {t) - Y k (s)J k (t)\t < Ke<^ [l + J-J , 

\Y k (s)J k+1 (t) - J k (s)Y k+1 (t)\ (l + 'pj ' t < Ke^ (l + 'pj ' . 



The first inequality follows directly from the second inequality of (41), and the second inequality is proved in 
the same way as the second inequality of (41), except that at the boundary where i = swe have 

\Y k {s)J k+1 (s) - J k (s)Y k+1 (s)\ = 2/ns. 

The details are omitted. 

It is also clear that the estimates above are uniform in A for A € (Ao/2, 2Ao). 

The estimates for Y k follow by the same estimates, since it is just a matter of multiplying each side of the 
inequalities by the factor (1 + k 2 ). □ 

We are ready to prove that there exist time-dependent exponential dichotomies on J + for the full system (5). 
First we define the spaces X s := X s x X s and y s := Y s x Y s . Note that y s C X C X s . As before, we can 
decompose those spaces into X s = ©feeZ^I an d y s — ©fcezJ^ with X k = Xf, x Xji and yf, = Y k s x Y k s . 
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For s > S3 = ra, we have that r = s and hence, for s £ [S3, 00), the system (5) reduces to 



U' = A(s;X)U with A(s;\) = 



/ 1 \ 

-\d 2 -i 1 

1 

V a -y 



(43) 



We consider the system (43) for s £ [j"i, 00) and record that (43) and (5) coincide on the smaller interval [S3, 00). 
The exponential dichotomy for (5) on the whole interval </+ will follow from the fact that the systems (5) 
and (43) are linked by the smooth transformation r(s) of the independent variable on the compact interval 
[si,S3]. Thus, since we are using r = s in (43), we see that if U(s) is a solution of (43) for s £ [ri,oo), then 
U(s) = diag(l, r'(s), l,r'(s)) U(r(s)) is a solution of (5) for all s <G [si,oo). Recall that there are constants 
< c < C such that c < r'(s) < C for all s £ R, thus dichotomy results for U will immediately give similar 
dichotomy results for U. 

It is easy to check (similarly to Lemma 3) that A(s; A) : X s — > X s is closed and densely defined with domain y s 
and that A(s; A) : y s — > y s is closed and densely defined with domain Z s x Z s . 

The Fourier coefficients of U = (ui, U2, W3, U4) T satisfy the system 

u[ = u 2 , 

k 2 1 

u' 2 = —Ul u 2 + u 3 , 

s s (44) 

u' 3 = M4, 

k 2 1 

u' 4 = Xui H — -u 3 u 4 , 

s z s 

of ODEs, where we omit the subscript k. We denote by $fc(s,<) the evolution operator corresponding to the 
system (44) and consider this evolution operator in either X s or 3^ s - 

Now we can use the earlier dichotomy results to show the existence of a uniform (s-dependent) exponential 
dichotomy for the system (44) and hence for (43). 

Lemma 14. There exists an eo > such that for any e £ (0, eo) there exists a K > such that for any 
A £ (Ao/2, 2Ao) there exists an s-dependent exponential dichotomy of (43) on J + such that the evolution operators 
$^(s,t;A) and $™(s, t; A) solve (43) and 

lin(*.*^)lk^) = ll*+(*.*;^)IUo'*j") < ^fH Al/4 -^- 4 \ s>t>n, 

ll$ c +"(M; A )ll£(;^) = ||*?'(s,t;A)|U C yt,y.) <Ke< l - s \ t > s > n. 

T/ie dichotomy is smooth in A /or A near Aq. 



Proof. Since Af s = and each is mapped into Xg under the flow of (43), we write for U £ X 1 

<$>%{s,t)U := ^$£(s,t)LV fc \ s > t > n, 

fceZ (46) 
$™(s,i)[/ := ^ $ c fe u (s,t)C/ fe e 2fc -, t > s > n. 
fcez 

Moreover, for each fc £ Z the evolution operator $fe(s, t) associated with (44) can be expressed in terms of <fik(s, t) 
and ipk(s,t). Indeed, it can be seen that 

Similarly, 

2 l/a#(m) ^(*,*) y 
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and 



2 \yA($(M)-^(M)) 4>%(s,t) + M*,t) , 



Introducing the temporary notation U k = (u,v) £ Xj,, where u{u\,U2) an d v = (1*3,1*4) , we have 

2 



wn(s,t)Uk\\%s = ](i+x) 



(47) 



Since we will take the supremum over all U k € Xj. such that ||f/fc||#i = 1, we may without loss of generality 
assume that v = \/\u, since all other choices will result in a smaller value of the right hand side of (47). For 
any such u and v € Xj., the condition that |j£4||^ t = 1 implies that ||w||l~t (1 + A) = 1. We therefore have 



sup \\3>%{s,t)U k \\ 2 



1 + A 



= sup — 

ll«ll^ t =l/(l+A) 

x k 

= 1 SU P Il0fe(s,*)u| 



X: 



which shows that 



Likewise, 



where 



As above, 



Il*fc(*.*)ll/:(^ 1 *i)=^ll0fc(*»*)ll/:(i*,jc i )- 
^ x l x i) + 



-y/\ipk(s,i) ip k (s,t) 



iM«,t) \ 

V* and = *£*(«,*)-*£(«,*)• 



ll**(*,*)IU(^)=^ll^( fl ' i )ll£(i 



\\n(s,t)\\c(xt,X£)=^Hk 



From Lemmas 12 and 13 it follows that 



ll<^M)ll£(^) ^e-i^M, 8 >*>n, 
\\*?(°>t)hixix S ) <Ke^- s \ t>s>n. 



By (46) and (48) we see that for s > t > r x 

W+M\\c{ X *, X °)= sup W{ Sl t)U\\ 2 Xs = sup \\j2n(s,t)Ukt 

\\u\\ xt =i \\n x t=v [ t£ 



i k ■ 



sup 5^|$jE( a ,t)E/ fc < sup 



1. ^ SU P u $ fc(M)IL 

** 11^11^=1 fcez 

2 p -2(A 1 / 4 - e )(s-t)||f), ||2 ^ = i ^2 e -2(A 1 / 4 -e)( S -t) 



ll^ll 



•V, f 



< sup V^e-^"-'^-^!!^! 

\M x t=l keZ 



(48) 



A similar calculation shows that $™ satisfies the second equation of (45). The estimates for Y£ follow by the 
same estimates since it is just a matter of multiplying each side of the inequalities by a factor (1 + k 2 ) 2 . 

Finally, the smoothness in A follows from the implicit function theorem [4, Corollary 3.1.11]. As A(s; A) depends 
linearly on A, we get that <j>^ cu (s, i; A) satisfies 



^< C "(M;A) = L4( S ;A ) + (A- A )B ] < cil ( 



where Bq is defined by (19). 
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For A € (Aq/2, 2Aq) and sufficiently close to Aq, the pair ($+(s, t; A), $™(s, t; A)) satisfies the fixed point equation 



$i(s, t; A) = $i(s, t; A )+$;(s, t; Ao)*+ t; A) + (A - A ) 



- f <f s + (s,T;X Q )B Q ^(T,t;\)dT 
Jn 



/s poo 
& + (s, t; Aq) B $ s + (t, t; X)dr- J ^(s, r; A ) B $ s + (t, t; A) dr 

s>t> n; 

$™(s, t; A) = $; u ( S , t; A )-$;(s, t; A )$ c + u ( S i , t; A) + (A - A ) f $' + (s, r; A ) S $?(r, t; A) dr 

.Jri 

/t poo 
$f(s,T;\ )B $™(T,t;\)dr + J $™(s,r; A ) B *+(t,«; A) dr , 

t > s > n. 

This fixed point equation is considered as a mapping on X s x X cu , where X s ' cu are defined as 

X s = U s {s,t)£L{X\X s ); W{s,t)\\x*= ™p e« x °M 1/4 -^ s -V\\$ s (s,t)U(x*,x*) < 

{ s>t>rt 

X cu = U cu {s,t)^C{X\X s )- ||$ C «(M)||^ = sup e- 2e (*-)||* -(* J t)|| jC( ^ i ^ ) <oo). 

L t>s>r t J 

With the estimates derived before on $i(s, t;A) and $™(s,i;A), it is easy to check that the right-hand side 
is well-defined in this space. With the implicit function theorem, it follows immediately that the mapping is 
smooth for A near Aq. □ 



6 The unperturbed adjoint system 

The dual space of X is X' = H~ 2 x H~ x x H^ 1 x L 2 (using the L 2 pairing). For the space X' , we make the 
decomposition 

X = ®k&X' k , 

where X' k are 4-dimensional pairwise orthogonal subspaces span{((a, b, c, d)e lk ') ; a,b,c 7 d € C} C X' . For 
W E X' k and U € Xk, we have the pairing 



(W, U) := WiU\ + w 2 U2 + W3U3 + W4U4, 



(49) 



where the bar denotes the complex conjugate. This means that we may use the standard inner product on C 
when computing the adjoint equation. 

Similarly, the dual space of X := H 1 x L 2 x H 1 x L 2 is X' = H^ 1 x L 2 x H~ x x L 2 (using the L 2 pairing) and 
we can make the same decomposition as above, i.e., 



X' 



1k&X' k , 



where X' k are the same 4-dimensional subspaces as above but are now regarded as subspaces of X' . For W £ X' k 
and U € Xk, the pairing is again as in (49). 

At the end of Section 4, we have investigated the solutions of the unperturbed linear system U' = A(s; Ao, 0)U 
on J_. The adjoint unperturbed system is 



W = ~A{s; A , 0)*W = -(At + B{s; A , 0)*) W. 



(50) 



Just as in the case of the unperturbed linear system itself, expanding W in a Fourier series shows that the Fourier 
spaces X' k are invariant under the flow of the adjoint system (50), and the Fourier coefficients satisfy the adjoint 
equation of (14), i.e., 

W' k (s) = - \A-(k)* +B(s;\,0)*] W k (s). (51) 
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It is well known and straightforward to check that the pairing of a solution of a linear system with a solution 
of its adjoint is constant. For our systems, this means that any two solutions U k (s) of (14) and W k (s) of (51) 
satisfy 

£(W k (a), U k (s)) = 0, and thus (W k (s), U k {s)) = (W k ( Sl ), U k ( Sl )) for any s G R. (52) 
as 

From [17, p. 56] it follows that if a finite-dimensional linear system has an exponential dichotomy on an interval J 
with constants K, n s and k u , then the adjoint system has an exponential dichotomy on J with the dichotomy 
constants K, — k u and — k s . Furthermore, if we denote evolution operators corresponding to the exponential 
dichotomy of the adjoint system (51) by $ s k± (s,t) and & k (s,t), respectively, then <J?|(s,t) = & k (t, s)* for t < s 
with s,t G J and %(s, t) = $ s k (t, s)* for s < t with s, t G J. 

On J_, the dichotomy constant K in (17) is independent of k, and so we immediately get the following estimates 
about the solutions of the adjoint system in the Fourier spaces X' k with norm 

WW II 2 — IIWp*II 2 I^i| 2 Kl 2 |w 3 | 2 , ,2 , 

Lemma 15. There exists a K > swc/i £/ia£ /or every k G Z\{0} and G 

||*jE(«,*)Wfc||jci < if e^l^ HWfcHxi, t<s< Sl , 
\\ms,t)W k \\ xi < Ke^ s ~^ \\W k \\ X ' k , s<t< si. 

Furthermore, for any solution W k (s) with W k (s\) G Ran($^_(sx, Si)), we have 

\w k (s)\ < K e W s -^\\W k ( Sl )\\x> 

for all s < si, where w k (s) denotes the fourth component ofW k (s). 

Similar arguments give the dichotomy of the adjoint system on J + . From (48), it follows that the solutions 
of the linearised system in the Fourier spaces 3>|(s, t) and $™(s,i) have an exponential dichotomy with a 
uniform constant K. The dual space of X s is denoted by (X s )', and for s fixed, this space is equivalent to 
x L 2 x H~ x x L 2 . The dual Fourier space is denoted by (X k )' ', and its norm is 

\\W k \\f x , y := \\W k ^\\\ xay = -^Kl 2 + Kl 2 + p^KI 2 + Kl 2 . 
On J_|_, we have the following estimates: 

Lemma 16. For every e > 0, there exists a K > such that for every k G Z\{0} and W k G X' k 

\\n s {s,t)W k \\ [X sy < Ke< s ~^ \\W k \\ (xty , Sl <t< s, 

\\ms,t)W k \\ (X sy < KeW-'X'-Q \\W k \\ {xty , Sl <s<t. 

Moreover, for any solution W k (s) with W k (s\) G Ran($^. s (si, s±)), we have 

\w k (s)\<Ke< s -^\\W k ( Sl )\\ {xsiy 

for all s > si, where W k (s) denotes the fourth component ofW k (s). 

Next we look at the adjoint system associated with k = 0. In Lemma 10, we have seen that the solutions space of 
the linear system at k = is spanned by Uqj, j = 1, . . . , 4. Now let Zq^ (s) be solutions of the adjoint system (51) 
with k = such that {•^o,j( s i)}(=i,- -,4 i s a dual basis of {?7o,j( s i)}j=i,...,4 (i- e -j {Zo,i{s\), Uqj(si)) = Sij). 
With (52), this implies 5ij = (Zq^(s),Uoj(s)) for any s < s\. The unbounded solutions C/0,2 and f/o,4 arc 
not unique, but they can be chosen such that Zq^ and Zq^ are bounded on J_, whereas Zq^i and Zq^ grow 
algebraically as s — > —00. A convenient choice for Z$^ and ^0,4 is 

ZoA*) : = i = 2 ' 4 ' (53) 
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where U = (— 114, U3, — M2, u{) if U = (ui, it2, W3, U4) 71 . It is easy to check that Z^(s) are solutions of the 
adjoint system (50) using that Uqj-i(s) are solutions of the original system (5). 

As will be shown below, a similar exponential dichotomy on J_ as in Lemma 15 also holds with norms in X' . 
Lemma 17. There exists a K > such that, for every k £ 7\{0} and W k € (X k )' , we have 

\\& k (s,t)W k \\ x , < Ife-Wt 8 "' 1 \\W k \\ X ', t<s< Sl , 
\\$t(s,t)W k \\ x , < Ke^ s ~^ \\Wk\\ x ,, s<t< Sl . 

For any solution Wk{s) with W k {s\) £ Ran($^(si, Si)) and wii/i Wk{s) denoting the fourth component ofW k {s), 
we have \w k (s)\ < Ke^ s - a ^\\W k (si)\\x' for all s < Sl . 

Proof. First we will show that the linear system (14) has an exponential dichotomy in X k . The proof is very 
similar to the one in section 4 with a slightly modified matrix M k ■ Define the matrix M k whose columns consist 
of eigenvectors of A_(k) that are scaled different to those in M k : 



M k = 



f-l/\k\ l/\k\ \_ 

1 10 

-l/\k\ l/\k 

V 1 1 / 



As M k consists of eigenvectors of A_(k), it follows immediately that A_(k) = M k D k M k . It is also straightfor- 
ward to verify that M k is a homeomorphism between C 4 and X k with ||Mfe||£(c 4 ,Ar fc ) — ^ \/2 as |fc| — > 00. 

Using the same ideas as in the proof of (17), with M k replaced by M k and exploiting the observation that 

\M^B(T/\k\;\ ,0)M k \ < f^!- Bup{i,|A-flf(a)|}, 

we find that 

\\n(s,t)U k \\ x <Ke-^ s -^\\U k \\ x , t<s<8i, 
M(s,t)U k \\ x < Ke^-V \\U k \\ x , s<t< Sl 

for any U k € X k for some constant K that is independent of k. As X k is finite-dimensional, this immediately 
implies the estimates of the Lemma. □ 

Finally we will show that, for large values of k, the solutions of (51) are close to the solutions of the asymptotic 
system W' k {s) = — A-(k)* W k (s). Recall that we denote the spectral projection onto the eigenspace of A_(fc) 
associated with the positive eigenvalue |fc| by P k and the complementary projection onto the eigenspace of A -(h) 
associated with the negative eigenvalue — |fc| by P|. 

Lemma 18. For every e > 0, there exists an N £ N and a S > such that, for every \k\ > N, we have 

||%(s, Sl ) - el^'^ni^) < eeW s - s >\ for all Sl - 5 < s < s u 
Sl ) - e) k ^\P°y\\ c(XD < ee l fe K s ^), for all s < 8l . 

Thus, for \k\ > N and W k {s\) e Ran($^(si, s\)), we have 

\\W k ( Sl ) - (P^*W k ( Sl )\\ XL < e\\W k ( Sl )\\ x ,, 

\\w k ( Sl ) - (PnywMWw < e\\W k ( Sl )\\ xi . 

Proof. The evolution operator $^(s,si) satisfies 

n( S , Sl ) = e^ s -^(P£y - f e -l fe l( s - t )(P,")*S(i;A o ,0)*$^i, Sl )di 
+ r e^-^iP^Bit; \ o ,0ym, *i) 



(55) 
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From its definition (12), we immediately see that \\B(t; A , 0)*|| £ (^) < Ce 2t and \\B(t; Ao,0)*|| £ (^) < e 2t 
for some constant C, independent of k. Hence, with the dichotomy estimates from Lemma 15, we get that for 
s < si and k e %\{0} 



|%(s,si)-e |fc|(s " Sl) (Pfc)*||£(^) <CK / e- |fc|(s - 4) e 2t e |fc|( *- Sl) d< 

J — oo 

+ CK /"VK'-VeWe-'^di 



2 + 

It is now easy to see that we can choose > large enough and 8 > small enough such that the first inequality 
of (55) is satisfied. 

With the dichotomy estimates from Lemma 17, we get 

\m(s,si)-e ms - Sl) (P£y\\c(x L ) < ~^=r K f e -|fc|(-t) e at e |fc|(t-. 1 ) dt 

V A/ I 1 J — OO 



< CK e ifci(«-«i) ( c : i C 2 S1 



2 v / fc 2_ +T^ V(i + |*|) 

< OAe e |fc|(«-«i) 

It follows that also the second inequality of (55) is valid when TV is sufficiently large. □ 



7 Matching the core and far field solutions 

In the next lemma we show that u is an eigenfunction that belongs to an embedded eigenvalue A of C + p if and 
only if U is a solution of (5) such that U(si) <G RanP£(si; A, p) n (RanP^(si; A, p) ® Pl b (si;X, p)). 

Lemma 19. Let u be an L 2 solution of (2). Then the corresponding solution U(s) of the system (5) is bounded 
in X as s — > — oo and decays exponentially with rate A 1 / 4 — e as s — > +oo in the sense that for any e G (0, A 1 / 4 ) 
there exists a constant K > such that 

\\U(s)\\ x . < Ke-^ 1 '"-^ (56) 

for every s > s±. Conversely, if U is a weak solution of (5) such that \\U(s)\\x is bounded as s£ J_ and such 
that \\U(s)\\x' decays exponentially as s — > +oo (with any decay rate), then it corresponds to an H 4 solution u 
of (2). 

Proof. If u is an eigenfunction of (5), it belongs to H 4 (M. 2 ). By Lemma 2, the associated solution U(s) of (5) 
is bounded in X as s e J_. In J + , for s > S3, the system (5) reduces to (43) and the decaying solutions of this 
system are series formed by Bessel functions Kf-, Jk and Yfc. The decay of Jk and Yk is asymptotic to l/\/s as 
s —> +00, and so these solutions do not give rise to L 2 solutions of (2). It now follows from Lemma 14 that U(s) 
decays exponentially as s — > +00, in the sense that (56) holds. 

Assume that U is a bounded weak solution of (5) which decays exponentially as s e J + . We need to show 
that the first component of U which we denote by u belongs to H 4 (M. 2 ) when regarded as a function of the two 
variables (r, ip) in radial coordinates. As U is a weak solution of (5), by Definition 1, U € Lf oc (J; Y) HiJ 1 oc ( J; X). 
Also, ||E/||x is bounded on J_, and and hence U € L^ C (M_; X). From Lemma 2 we know that u g i/ 4 oc (K 2 ), 
so we only need to worry about the decay properties of u as r — ^ 00 (i.e. as s — »• 00). From Lemma 14 and the 
definition of Af s it follows that for any < e < eo and s > si 

||«WI| ia( 5i) <Ke-^ /i -^ s -^\\U(s 1 )\\ X s 1 
and so it is clear that u € L 2 (R 2 ). From (2) it then follows that u £ H 4 (M. 2 ). □ 
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Recall that is the radially symmetric eigenfunction associated with the embedded eigenvalue Ao when p = 0. 
Let Z7* be the associated solution of (5) with p = and A = Ao, defined for s € R, i.e. £/* = (it*, u*, Aw*, (Ait*)') T . 
Define X := H 1 x L 2 x H 1 x L 2 and recall that X = H 2 x H 1 x H 1 x L 2 . Let 

: = {C/G AT; P*(si; A 0l 0)U = U}, 
El: = {UeX; P^( Sl ; A , 0)U = U}, 
E c } : = span{C/o,i(si), C/ , 3 (si)} C X, 

where Uoj are defined in Lemma 10. Roughly speaking, El and El consist of the initial values at s± of solutions 
of (5) with p — and A = Ao which decay exponentially as s — > oo and as s — oo, respectively, and El © -Ef 6 
consists of the bounded solutions on J_. Note that the norm of X is used for while the norm of X is used 
for El and E c }. We have £^ n (El ffi Ef 6 ) = span{[/*(si)} since Ao is an eigenvalue of C with multiplicity 1. 

Next, we introduce a new Hilbert space X such that X C X C X , and special solutions Vk,j, k G Z, j = 1, . . . , 4, 
such that {Vfcj (si)} is a basis for AT. 

We have seen that the unperturbed system (5) decouples when p = so that the subspaces Xk and Xk are 
invariant under the flow of (5) with p = 0. For fc = 0, we pick Vb,i(si) := f/*(si) € E c } n Note that 
there are no other solutions in Xq which decay exponentially as s — > oo. We pick Vo,2(si) G E^ such that 
Vb,2(si) i E s + and E c } = span{Vb,i(si), F , 2 (si)} (thus span{F ,i(si), V 0t2 (si)} = span{[/ ,i( s i)> ^0,3(si)})- 
Next, we will choose Vq^(si) and Vb, 4 (si) € X such that they belong to the span of £/o,2(si) and ^ r o,4(s2)-In 
order to do this, we introduce a dual basis Wo,j(si) of Vbj(si) and choose Wo^(si) := U^~(si), where U^~(si) = 
(— U4,(si), U3(si), — M2(si), Ui(si)) and Mj(si) are the components of t/*(si), while Wo,4(si) is any other vector 
so that span{ W , 3 (si), Wo,4,(si)} = span{Z , 2 (si), Z 0i4 {si)}. The remaining vectors Wb,i(si), W Q ^(si), Vb,3(si) 
and Vqa(si) are determined by the conditions that {Wo,j(si); j = 1, . . . , 4} and {Vcy(si); j = 1, . . . , 4} are dual 
bases: 

(W 0l i(«i),V 0l i(si)> = fy, 1 = 1,... A- 



We use the notation E™ := span{Vo,3(si), Vb j4 (si)}. We normalise the vectors such that ||Vbj(si)||x = 1 for 
j = 1, ... ,4 and note that, for V G Xq, we have H^Hx = ||V||;t:. Define Wqj(s) so that it satisfies the adjoint 
system (50) (and hence (51) with k = 0) and passes through Woj(s±) for s = s\. From (52) and the relation 
above, it follows immediately that (Wo,j («), Vq,i(s)) = Sji for all s < s\. Furthermore, from (53), we conclude 
that Wo,3(s) and Wo,4(s) are bounded solutions of the adjoint system on J_. 

Next we consider k ^ 0. The spaces Xk and Xk are four-dimensional, and Si n Af^ and E+ n Afc are one- 
dimensional, El n X fc and £™ n X k arc two-dimensional, and ElC\E s + C\ X k = {0} = £! n £^ n 4, since the 
multiplicity of the eigenfunction [/* is 1. Using this, we define base vectors in Xk and Xk as follows: For k ^ 
we pick Vfc,i(si) G (and hence Vk,i{ s i) ^ E"). We also pick Vk^isi) and V^^Si) so that they belong to El 
(and hence do not belong to E!L). Thus {Vk,i{s\), Vk^{si), Vk^(si)} span a three-dimensional subspace in the 
four-dimensional spaces Xk and Xk- We normalize the solutions Vk,j such that for k G Z\{0}: || Vfe^si)!^ = 1 
and ||Vfej(si)||x — 1 f° r j = 2,3. Hence there exists a unique (up to multiplication by a unimodular constant) 
vector Wk,4 G such that 

(W fe ,4 , V k ,j) = for j = 1,2,3 and ||W M (si)||^ = 1. 

Then (Wk, 4 (si), V) = for V G Ran($| + (si, Si))+Ran($£ _ (si, Si)) and hence W M (si) G Ran($£" + (si, s x )*)n 
Ran($)* _(si, si)*) = Ran($£ s + (si, si)) n Ran($'^ _(si, si)). We take the one remaining solution in Xk and Xk 
such that (W M (si), Vi,4(*i)) = 1 and ||T4, 4 (.si)|U = 1. Then F M (si) El U E s + as (W M (si),^) = for 
V e El + E%. 

Let X be defined by 



X := I U = 



^ fed 7. fcf=K * kPW. X ) 



j=i,-A 



fee 



i=2,3 
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Note that X is the direct sum of two Hilbert spaces, which are closed subspaces of X and X, respectively. It 
follows that X is a Hilbert space. 

Note that E^_ and El are both closed subspaces of X: Indeed, 

E\ = clyspan{F fej i(si)} fc6 z = cl x span{14 j i(si)} fceZ 

El = cl x span{F fcj2 (si),Ffc,3(si)} fceZ \ { o} = clx span{F fc , 2 (si), V ki3 (si)} keZ \ {0} . 
It is clear that is a closed subspace of X since it is finite-dimensional. 
Define i : E s + x {El ® E°}) xlxK^Iby 

i{U^ ttf + C/ ch ; A, p) : = P"(si; A, p)C/ " + P! 6 ( Sl ; A, p)C/ cb - P+(si! A, flCfl 
= P"(si; A, p)Etf + P! 6 ( Sl ; A, p)C/ cb - P+(si; A, 0)[/ o s , 

where wc recall the definition (30) and note that the last equality holds since p(s) = for s > s\. We will see 
that the range of l is a subspace of X, and that l is smooth into this space. For this we need a more explicit 
formula for l. By (30), (31) and (49) evaluated at s = t = Si, we have (see (19) for the definition of Bq) 

Pl b ( Sl ;X,p)U^ = [/ C& + r ^( Sl ,T;X ,Oy(T) 2 (X-X a -p(T))B U c _ b (T;X,p,U^)dT (58) 

J — oo 

PH(si;X,p) = P^( Sl ;X 0l 0)+ f 1 ^( Sll T;X Q ,0y(T) 2 (X-Xo- p(r))B Q ^l(r, Sl ;X,p)dT (59) 

pOO 

Pl(s V ,X,0) = Pl( Sl ;X 0> Q)- <S>l u ( Sll T;X Q ,Oy(T) 2 (X-Xo)B <S>l(T, Sl ;X,0)dT, (60) 

J si 

where we recall that r(r) = e T for r < S2, so that r'(r) 2 = e 2r in this interval. Hence we may write 

/>oo 

L(U$,UZ + UZ b ;X,p) = K + UZ b -UZ+ / $™( Sl ,r;Ao,0)r'(r) 2 (A-Ao)i?on(^si;A,0)[/ s dT (61) 

J si 

$ c _ s ( Sl , r; A , 0)e 2T (A - A - p(r))B [U c _ b (T; A, p, C/ cb ) + *-(r, «i; A, flttf] dr. 
Lemma 20. The map l:E s + x(E^® E c }) xixK^I is smooth. 

Proof. We have seen in Theorem 2 that (Utf,X,p) i-> P"(si; A, /5)[/q and (U§ b ,X,p) i-> P™{ Sl ; X, p)U§ b are 
smooth as functions from x R x K to X C X and from F'xlx'KtoICl, respectively. 

Hence it suffices to prove that Pf (si; A, 0) is smooth from E+ x R to AT. We do this by studying the terms of 
(60) separately. It is clear that L^f G X. 

Next, we study the integral term, and note that by Lemma 14, for r > si, $1(t, Si;A,0) : # S1 — > X T with 
norm bounded by X e ( Al/4 - £ )( r - s i) and $^(si, t; A , 0) : J T with norm bounded by Ke £(r ~ Sl) . Recall 

that X s = H 1 x L 2 x H 1 x L 2 and CP = H 2 x H 1 x H 2 x H 1 with s-dependent norms. Thus B : X T -> ^ r 
is bounded with norm r. Using the exponential estimates of Lemma 14 and that y si C X C X C , there 
exists a constant C > such that for any e > sufficiently small 

/>oo />oo 

/ $™( Sl ,r;A ,0)r'(r) 2 S $;(r, Sl ;A,0)t/ s dT _< C / *™( Sl , r; A , 0)B o $ s + (t, Sl - A, 0)C/ S 

<CK 2 ^ xl/i - 2 ^ / re-^'-^drWU^ 

J si 

- (Ai/4_2e) 2l|C/ ° S|lx ' 

for some constant C. To show that the integral term is smooth in A into X, it suffices to prove that for n > 1 

y *o> 0)r'(T) 2 B — <Z> s + (T, Sl ; A, 0)[/ o s dr 

belongs to X. This follows since, by Lemma 14, -73-5-$+ (t, Si; A, 0) satisfies a similar exponential decay estimate 
as $^_(r, si; A,0) docs. It follows as above that the integral in question converges in y si C X. Smoothness in 
Uq is immediate, since l is bounded and linear in Uq into X. □ 
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Lemma 21. The operator A 2 + 6 + p has an embedded eigenvalue A > if and only if there exist Uq G E+, 
Utf G El and U§ b G EL b such that 

i(U^U^ + U^ b ;X,p) = 0. (62) 

Proof. If A is an eigenvalue of C + p, then by Lemma 19, the corresponding solution of the system (5) is bounded 
as s — > — oo and decays exponentially as s — > +oo. Hence there exists a solution of (5) with initial condition 

P+(si; A, p)Uq = P"(si; A, p)Uq + P cb (s i; X,p) 

at s = si, i.e. (62) holds. 

Conversely, suppose that (62) is satisfied for some (C/ S , lift + U§ b , A, p) G EL x (£^ ijf. 6 ) xlxi Then there 
exists a solution of (5) with initial condition Uq = Uq + U§ b . By Lemma 19, this implies that A is an eigenvalue 
of C + p. □ 

Lemma 22. The subspaces EL b ®EL and EL have complements in X denoted by EL a (BEL and E™. Moreover, 
(EL a © EL) H EL U is infinite- dimensional and has a basis with elements Vo,3(si), Voa(si) G Xq, 14,4 (si) G X k , 
k G Z\ {0}. 

Proof. Recall that E<L a = span{Vb, 3 (si), Vb, 4 (si)}. Let 



: =span{V r M («i),V r M (si); fc G Z \ {0}}, 
E? : = spajy{V k , 2 ( Sl ),V k , 3 ( Sl ),V kA (si); k G Z}, 
where the closures are taken in X. It is easy to see that these spaces have the desired properties. □ 

Let Q be the projection in X onto Rani(-, •; A , 0) = EL + (EL © EL b ) such that 

kcrQ = E^Ci{E m ® EL). 

Note that RanQ and kcrQ arc closed subspaces of X, and Q is therefore continuous. 
Equation (62) is equivalent to the pair of equations 

Qi(U^ + U^;X,p) = 0, 
(I-Q)l(U o s ,US + U o J} ;\,p) = 0. 

Lemma 23. For (A, p) in a neighbourhood of (Ao,0) GlxR, the first equation of (63) has a unique (up to 
constant multiples) nonzero solution (Uq,Uq + UQ b ) which depends smoothly on A and p in this neighbourhood. 
We write U§(X,p), U$(\,p) and U Q ab (X,p). Furthermore, U^(X o ,0) = C/»(si) = t/ o cfc (A o ,0) and Utf(\ ,0) = 0. 

Proof. For (A, p) fixed, Ql is a linear mapping from EL x (i?" © -Ef 6 ) to RanQ. It is clear that 

kerQt(-, ■; A ,0) = span{([/*(si), L7»(si))}. 

By Lemma 20 and since Ran Q is closed, it follows that Ql is a smooth mapping in its arguments. Let D be an 
affine hyperplane of EL x (EL © EL b ) such that L> n span{(£/*(si), 17, (si))} = {(*7*(si), *7*(si))}. The implicit 
function theorem then implies that for (A,p) close to (Ao,0) the first equation of (63) has a unique solution 
(U^U^ + U^) = (U s (X,p),U${X,p) + U§ b (X,p)) G D in a neighbourhood of (J7*(si), {7*(s x )). Moreover, C/J, J7 " 
and U§ b are smooth in their arguments. □ 

For (A, p) in the neighbourhood obtained in Lemma 23, we let 
F(X,p) : = (I- Q)l(U^(X,p),US(X,p) + U£ b (X,p)) 

= (I-Q) f 1 $-( Sl ,r;Ao,0) e 2r (A-Ao-p(r))i?o [^'(t; X, p) + $L(t, Sl ; X, p)U£(X, p)] dr 

•J — DO 
/>OG 

+ (I-Q) $™( Sl ,r;Ao,0)r'(r) 2 (A-Ao)i?o$;(T, Sl ;A,0)C/ o s (A, / 5)dT 
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where C/q & (s; A, p) corresponds to F/!: b (s, A, p, [/g b (A, p)) so that UQ b (s; Xq, 0) = C/*(s). We see that solving (63) is 
equivalent to solving F(A, p) = 0. 

On Ran(7 — Q) C X, the X-norm is the same as the A'-norm, and so we solve F(A, p) = in X. For k £ Z \ {0} 
let F k (X,p) := (iy fe , 4 (si),F(A,p)), and for fc = and j = 3,4 we let F 0tj (X,p) := (Wo,j(si), F(X, p)). Define 
Wk^{s) so that it satisfies the adjoint system (50) (and hence (51)) and passes through Wu,i{si) for s = si. As 
Wfc, 4 (si) £ Ran($^ + (si,si)*)nRan($| _(si,s 1 )*) = Ran($^ + (si, si)) D Ran(%_(si, si)), we get for k ^ 



F fc (A,p) = / (W M («i), $i s ( Sl ,r;Ao,0) e 2r (A-Ao-/5(r))B [F cb (r; A, p) + *«(r, sv, A, /5)C/ "(A, p)] ) dr 

*/ — OO 

/>oo 

+ / (W M (si), *^(ai,r; A , 0)r'(r) 2 (A - A )B $;(r, Si; A, 0)t/ o s (A, p)) dr 
(W kA (r), e 2T (A - A - p(r))B [EtfV; A, p) + *^(r, si; A, p)C/ «(A, p)] ) dr 

) 

/•OO 

+ / (W M (r),r'(r) 2 (A-Ao)S $;(T, Sl ;A,0)C/ s (A,p))dr 

J si 

and for k = and j = 3, 4 we have similarly 

F ,j(X,p) = [ S \woAT),e 2 ^X-X a -p(T))B [UZ b (T-X,p) + <Pl(T, Sl ;X,p)U^ (64) 

J — OO 

+ r ( W °,i ( r )' r ' ( T ) 2 ( A - Ao)S $; (r, Sl ; A, 0)C/ S (A, p)) 



Lemma 24. For (A, p) in a neighbourhood of (Ao,0) £ R x 7?., £/ie equation (63) /ia,s a nontrivial solution 
(Uq(X, p), F/^(A, p) + F cb (A, p), A, p) £ F* x (£« Ff 6 ) x R x K if and only if F k (X, p) = for k £ Z \ {0} and 
F 0ii (A,p) = 0/ O ri = 3,4. 

Proof. Suppose that (F S (A, p), t^(A, p) + F cb (A, p), A, p) £ F^ x (F« F^ b ) xlx^ solves (63). It is then 
clear from the definition of Ffe that F k (X, p) = for k £ Z \ {0} and Foj(A,p) = for j = 3,4. Conversely, let 
Ffe(A, p) = for £ Z \ {0} and Fo j(A, p) = for j = 3, 4. By Lemma 23, the first equation of (63) is satisfied, 
so it remains to check the second equation of (63). Recall that the basis vectors in Ran(7 — Q) arc Vkj(si), 
where j = 4 for k ^ and j = 3, 4 for k = 0. The coefficients of (I — Q)l(Uq, Uq + U^; A, p) with respect to this 
basis are then Fqj (j = 3,4) and F k , k £ Z\ {0}. Since all these coefficients vanish, the conclusion follows. □ 

Lemma 25. The equation Fq^(X, p) = defines X as a smooth function of p in a neighbourhood of p = such 
that A(0) = Ao. Furthermore, 

A , (0) _ = _ O^oW^ 
X!^ w»(r) 2 r(T)r / (r) dr 

where p is the Fourier coefficient of p corresponding to k = 0. 

Proof. By Lemma 20 it follows that F0.3 is a smooth function of A and p in a neighbourhood of (Aq,0). Note 
that 

ap roo poo 

l ' 3 (A ,0)=/ (Wo,3(r),r'(r) 2 Bo^(r))dT = / »o, 3 (t)«,(s)^W 2 d», 



<9A 

where we have used that F o s (A o ,0) = U*(s) = F cb (A ,0) and U$(X ,0) = 0. We recall that W 0)3 (si) = E^(si) 
and that Wo, 3 (s) satisfies the adjoint system (50) for s £ R. It can be verified that 

Wb l3 (s) = 

r{s) 

for s £ K, where U = (—1*4, W3, — U2, Wi), and Uj are the components of {/, j = 1, . . . , 4. Hence 

SFq.3 , 



dX 



(A o ,0)=/ u*(s) 2 r(r)r'(r)ds > 0. (65) 
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The last inequality follows since the integral is positive (using that u* is an eigenfunction). 

Since dFo t s/dX(Xo, 0) 7^ 0, we can solve the equation Fo^(X,p) = by the implicit function theorem for A as 
a function of p, and this solution is a smooth function X(p), defined in a neighbourhood of p = 0, such that 
A(0) = Aq, and A'(0) is given by 



A'(0)p = -^(A 0) 0)p/^(Ao,0) 



fllu^rfp^e^dr 



IZo U *( T ) 2r ( T ) rl ( T ) dT 

as claimed. □ 

Since we have solved i*o,3 = for A in terms of p, the remaining equation corresponding to k = is Fq A (X(p), p) = 
0. We define 

G (p):= r (Wo^e^A^-Ao-pCr))^ dr 

J —OO 

poo 

+ / (W 0A (r)y{T) 2 {X{p) - X Q )B ^ + (T,s i; X(p),0)U^X(p),p)) dr, 

J si 

and for k 7^ 0, 

G k (p): = F k (X(p),p) 

(W kA (r),e 2T (X(p) - A - p(r))B [[/ cb (r; X(p),p) + <f u _(r, Sl ;X(p), p)U% (X(p) , p)] ) dr. 

pOO 

+ / (W kA (T),r'(T) 2 (X(p)~X a )B ^(T, Sl ;X(p^,0)U^X(p),p))dr. 

J si 

Lemma 26. The mapping Q : 1Z — > If defined by 

G(p) = {G k {p)} k ei 

is smooth. 

Proof. We first verify that the range of Q belongs to l\. To do this, we split the expression for G k {p) (k 7^ 0) 
into three terms, which we deal with separately: 



G k (p) =(A(p) - A ) / (W M (r), e 2T U cu (r)) dr - / (W M (r), e 2T p(r)U cu (t)) dr 

J — OO J— OO (nn\ 

f oo (66) 
+ (X(p) - A ) / (W kA (r), r'(r) 2 U s (r)) dr, 



where we used the notation 

U cu (t) := [US b (r;X(p),p) + ^ u _(T, Sl ;X(p),p)U^X(p),p)}, 
U s (r) :=$l(T,s i; X(p),pM(X(p),p). 

Then B U cu (t) e {0} x {0} x {0} x H 2 ^ 1 ) and B U s (t) e {0} x {0} x {0} x H^S 1 ). Furthermore, by its 
construction, we have that W kA {s\) € Ran($^ _ (si, s\)) n Ran($£ s + (si, s\)) for all k £ Z\{0}. Thus Lemma 16 
implies that for any e > there exists a constant K such that for every k £ Z\{0} and s > s\, 

||W M (*)|| C *.)' < Ke^ s ~^ 

as the norms on X' and (X Sl )' are equivalent and || VFfc^si)!!^' = 1. 

Now observe that BqU s (t) vanishes for all components except the last one, so we only need an estimate on the 
last component of W kA (s), which we denote by w kA . Then the estimate above gives that there exists a constant 
K independent of k such that 



\w k 



:a(s)\ < Ke< s ' Sl \ for all k € Z\{0} and s > s u (67) 
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as X' = H- 1 x L 2 x H' 1 x L 2 . Similarly, from W kA {si) G Ran($£ _(si, si)), Lemma 15, X' = ff" 2 x ii"" 1 x 



-fc,- 

r-i ^ r2 



H x L and || Wfc ) 4(si)||x' £ II W^fc, 4(si)||^' = 1, it follows that there is some constant K such that 

kfc, 4 (s)| < ife |k|(s " Sl) for all fc e Z\{0} and s < Sl . 

First we look at the last integral in (66). Let w|(t) be the first component of the fc-th Fourier coefficient of 
<&+(t, si; X(p), p)Uq(X(p), p), then the definition of £?o in (19) gives 



(WmW/W'So^'W) dr= / r'(r) 2 u; fc ,4(r)u|(T)dr 



J Si 

From Lemma 14, it follows that, for any e > and r > si, 



\\$l(T, Sl ;\(p),p)UZ(\(p),p)\\ x . < Ke-^ yi -^-^\\UZ(\(p),p)\\*°i- 
This implies for the Fourier coefficients ui that 

E + |^(T)| 2 <^ 2 e- 2 ^ 1/4 - e )^--)||C/o s (A(p),p)||^. 



fcGZ\{0} 

Combining this with (67), we see that for any p, we have 



E (l + fc 2 )(r/(r) 2 U ; M (r)^(r)dry<C £ (1 + fe 2 ) ( f°° Ke<^\u%{r)\ dr) 

GZ\{0} ^ Sl ' /cGZ\{0} ^ Jsi ' 

^ C E CT e " e(T ~ SlW ) (/°°^--)(l + fc 2 )|^(r)| 2 dr' 



fcez\{o} v 31 7 x 31 

<-fe 3£(T - Sl) E 2r 2 (f + ^)|^(r)| 2 dr 

poo 



where C = C(e) denotes the different constants occuring. 

Next, we look at the first integral in (66). Let it™ be the first component of the fc-th Fourier coefficient of 
U§ b (r; X(p),p) + (t, si; A(/5), p)Uq (X(p) , p) . The definition of i?o gives that the first integral can be written as 



si psi 

2r r>_ TTCUf^W J„ / „2t' 



(W kA (T),e 2T B U cu (T)) dr = / e 2r w kA (r) u™{r) dr. 

3 J — OO 

As <f>" leads to solutions with an A-norm that is exponentially decaying at — oo and is bounded in the 
X-norm, there exists a constant K such that the Fourier coefficients m™( t ) satisfy 



E (f + fc 2 )|^T(r)| < (l + fc 2 ) 2 |^fM 

fc<EZ\{0} fcGZ\{0} 



2 

< if 2 



for all r < si. Together with the fact that (68) implies that |u>fc i4 (r)| < K for all r < si, this gives 

E (l + k 2 )(^J S \ 2T ^ k ~J^)^{T)dT^j <K 2 E (1 + fc 2 ) (y" Sl e 2si e 2 ^ Sl M^f(r)|dr^ 



fcez\{o} s 7 fcez\{o} 



< ^e 2si Yi ( r e2{T ~ si)dT ) ( r (! + fc2 ) e2(r " si) i"rwi 2 ^) 

fcez\{o} 7 ^"'- 00 7 

ff2 2si psi jy-i 2s x 

-^^JJ E (i + ^ 2 )K"WI 2 ^<^^- 



feGZ\{0} 

Finally, let i/(r) be the first component of p{r)U cu (r), so that v(r) = p(r)u cu (r). As u cu G i? 2 (S' 1 ), its 
if 2 norm is uniformly bounded on (— oo, Si] and p G L 2 (J^\H 1 / 2 (S 1 ),e 2s ds), Lemma 6 implies that v G 
L 2 (J_; i/ 1 / 2 (S' 1 ), e 2s ds). 
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Denote the Fourier coefficients of v by v k . Then the second integral of (66) can be written as 



(W kA (T),e 2r p(T)U cu (T)) dr = / e 2T w kA (r)9 k (T)dT 

J — OO 

and the estimate (68) on the decay of Wk^ implies that there exists a constant C snch that 
e^^M^^drX < ( r K 2 e 2 ^ +1 ^- S A ( T e 2 ^(r) 2 dr 



\J — oo / w— oo 



oo 



and so 



Y, (! + fc2 ) if 1 e 2T w kA (r)^{r)dr\ <C £ (I + k 2 ) 1 ' 2 T e 2T d k {r) 2 dr 

fceZ\{0} V-oo fceZ\{0} 

1 2 



= C\W\\ L 2(j_. H i/2( S i^ e 2s ds - ) < OO. 

Hence the second term also belongs to I 2 , and so the proof of the claim that the range of Q is contained in If is 
complete. 

Smoothness follows since the integrands arc smooth in p and since the derivatives of arbitrary order of the 
evolution operators s±; X(p), p) and $"(t, s±; \(p), p) belong to the same exponentially weighted space as 

the evolution operators themselves (see Theorem 2 and Lemma 14). □ 

Finally we consider G'(0). Since U* is radially symmetric (and hence belongs to Xq) we have for k ^ that 

/•OO 

G' k (0)p= / (W k Ar),r'(T) 2 (X'(0)p- p(t))B U4t)) dr 



w kA (T)p k (T)u*(T)e 2T dr, 

> 

where p k is the fc-th Fourier coefficient of p. For k = we have 



S-oo w 0A (T)u*(T)r'(T) 2 dr f s ^ ^ 



G'o(Q)P=- w 0A {s) p (s)u*(s)e 2s ds ^ — — / p (s)u4s) 2 e 2s ds. 

J-00 J-oo U *( r ) r ( T ) r ( T ) dT 

To rewrite the preceding expressions, we define 



r] k (s) := Wfe,4(s)u*(s)x(_ 00 , Sl j(s) (69) 

for k £ Z \ {0}, and set 



?7o(s) := 

Then we may write 



7— r . , /!L w o,4 {t)u*(tY (t) 2 dr 

W AS) + MS) —Too M2 M <M , 

LooMt) r{Ty (t) dr 



w*(s)x(_ 00)S1 ](s). 



si 



g'(0)p=i- e 2T m (r)p k (r)dr ] 

I J-00 ) k! zi 

For any k <G Z, we have J] k e lk - e ft: indeed, (68) shows that |w m (t)| < Ke^-^ for any r < s x and fc € Z\{0}, 
while I 3 1 an d |u>o 4I are bounded on J_, so that there exists a constant C such that 



C 1 \Vk{r)\ 2 e 2T dr < sup K(r)| 2 C e (2|fe|+2)(T ~ Sl) dr < sup u.(t) 

J —00 r£{ — oo,si) J — 00 t£( — oc,si) 



2 ° , (70) 

2|fc|+2 W ; 

From the definition of £/'(0)/5, it follows immediately that G'(0)p = if and only if 

1 e 2r Vk (rfp k (r)dr = (71) 



00 



for all k € Z. Thus, if we define 



:= span{?? fe e ifc ^; fc S Z}, 

where the closure is taken in TZ, then it can be seen that M. is the orthogonal complement in 1Z of kcr Q'(0), and 
so it = ker^'(O)0M. 

37 



Lemma 27. Q'(0) : M. — > If is a linear homeomorphism. Furthermore, the spaces ker^'(O) and M are both 
infinite-dimensional. 

Proof. It is clear that G'(0) : R — > If is bounded since by Lemma 26, Q is smooth in a neighbourhood of 0. 
We need to investigate the subspace AL Let -q £ Ai be arbitrary, then 

» 7 (s,¥0 = X> fc »fr(s)e < **\ (72) 

The upper bound estimate (70) implies that 

= + fc 2 )i/> fc |2 f 1 h fe (r)| 2 e 2 ^r < C £ 



(73) 



Next we derive a lower bound for (Ml?,- Since u*(si) 7^ 0, there exist ?and S > such that ^(s) 2 > e 2 for every 
s £ (si — (5, si). Lemma 18 shows that, for fc large, W kA is close to solutions of the system at infinity, both in the 
X and X norms. This allows us to get a lower bound on ^^^(s)! for k large. Let e > and K as in Lemma 18. 
As W kA (s\) £ Ran(<&£(si,si)), it follows that ^{s, s\)W kA {s\) = W kA (s), and hence 

\\W kA (s) - e^ s -^(P^*W kA ( Sl )\\ xl < ee^ s -^\\W kA (si)\U' = ee^ s ~^ 

for |fc| > K . Thus we get for the fourth component w kA {s) 

\w kA {s)\ > e l fc l^ Sl )|((^)*^ M ( Sl )) 4 | - \w k (s) - e^-^((P^rW kA ( Sl )) 4 \ 

> e ^ s -^\((P°yW kA ( Sl )) 4 \ - \\W kA (s) - e^ s -^(P^rW kA ( Sl )\\ xl 

> e |fc|(-.0 (|((P«)*W M ( Sl )) 4 | - e ||W M (*i)|U0 > e l fc K s - Sl) (K, 4 (si)| - 2e||W M (.s 1 )|U0 
>el fc K s - Sl '(| U ; M ( Sl )|-2e). 

To get a lower bound on w kA (si), first note that Lemma 18 implies that ||Wfc^(si) — {P k )*W kA {s\)\\x' < e 
and that Ran(P^)* = span{(-|fc|, 1, 0, 0) T , (0, 0, -|fc|, 1) T }. Thus there exist a k , £ C and W k £ X' k with 
II Wit H*' < 1 such that 

W kA ( Sl ) = a k (-\k\, 1, 0, 0) T e* + /3 fc (0, 0, -|fc|, l) T e* + e^ fc . 
As Vfe.i(si) <S £4.) its first component has to be a multiple of K^(Xq r(si)) and hence 



V kA ( Sl ) = C k e 



ik- 



( K lkl (xl /4 ri ) \ 
r'( Sl )Xl /A K' lkl (\l /4 ri ) 

xl /2 K lkl (xl /4 n) 

\r'( Sl )Xl /4 K' lkl (Xl /4 ri )J 



C k e 



ik- 



x l/2 y I Withw fc 



^| fc |(Ao 7 Vi) 



r'( Sl )Aj /4 iT'(A^ 4 r 1 ) / 



where n = r(si) and Cfc is such that 1 1 T//c, 1 ( s 1 ) II -^f = L i.e. 



1 



(1 + Ao) 



(ft* + 1) (^(A^n)) + ^ (r'( Sl )^ fe| (Aj /4 rx)) 
Since (W kA (si), Vfe,i(si)) = 0, we have 

= C fc (-|fc|if| fe |(Aj /4 rx) + X 1 /A r'(s 1 )K' lkl (X 1 /4 r 1 )j (a k + f3 k ^X~ a ) + e(W k , V kA ( Sl )). 



(74) 



From (9.6.23) in [1], we see that K\ k \{z) > for any z > and (9.6.26) implies KLAz) = -if| fe |_i(z) - ^ K\ k \{z) 



for any z > 0, hence K! k , (z) < for any z > 0. So we see that -|fc|A'| fe |(Ao /4 n) < Oand X^r> (s^K^X^n) < 
0. A short analysis gives that 

-2/v/TTa^ < C k HklK^X^n) + Aj /4 r'( Sl )^|' fc | (A^n)) < - W2(l + A ). 
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Furthermore, Vfc^si))! < || Wk\\ X'\\Vk,i(si)\\x < 1, and we can conclude from (74) that a* = — /3fc\/Ao + 

0(e). Finally, ||W M (*i)|U' = 1 gives ^ (a 2 + ft) = 1 - 0(e 2 ), and hence (1 + A )/3 2 = ^ - O(e). Thus 
there exists a C > such that |iffc,4(si)| > C for all |fc| > K. This implies that there exists a C > such that 
\w kA {s)\ > Ce ms ~ Sl) for every s < s x and |fc| > 

Combining the lower bounds on u*(s) and 1/^,4 (s), we find that there exists a 5 > such that for |fc| > if 



? 2 C 2 



77 fc (r) 2 e 2T dr > ? 2 C 2 / e (2|*l+2)(r-0 dr = ( 1 _ e -(2|fc|+2)^ 

7 S1 _ 5 2\k\ + 2 \ 



2\k\ 



c 



> _l_r_(l_ e -M) > 



c 



(75) 



rsi 

for some positive ^-independent constant C. Since for all k £ Z, / rjk(T) 2 e 2T dr > 0, the constant C above 

«/ — oo 

%(r) 2 e 2r dr > 



can be modified so that 



also for \k\ < if. Hence it follows that 



(l + fc 2 )l/2 

| = E( 1 + fc2 ) 1/2 M 2 f ^(-) 2 e2T ^ > c Ei 



a fe | 2 . 



The upper and lower bounds on show that rj G *M if and only if 77 is given by (72) and {a/cjfcez G Z 2 . 

As we have seen that the mapping G'(0) is bounded above, it is sufficient to show that it is bounded below to 
conclude that G'(0) is a linear homeomorphism from M. to I 2 . From its definition, it follows that 



si 

2„2t 



g'(0)r,= ^a k J k(r)| 2 e 2T dr 
and so by (75) and (73), we see that 

l|0'(OMI?» =£(1 + *>*| 2 (T \^)\ 2 e 2s ds) 2 > C^Wkl 2 > C'\\r,\\l. 

It remains to show that the spaces kerC?'(0) and M. have infinite dimension. For A4, this follows directly from 
its definition. Next, consider the characterization of ker(J'(0) given in (71). We proved above that the functions 
Wk,i{s) that appear in the definition (69) of rjk satisfy |tUfc,4(si)| > C uniformly in \k\ > K 1 which implies that 
the space ker Q' (0) has infinite dimension as claimed. □ 

We are now ready to complete the proof of Theorem 1. 

Proof of Theorem 1. By Lemma 21, if (A,/5) is sufficiently close to (Ao,0) then A is an embedded eigenvalue for 
A 2 + 6 + p if and only if (62) holds. We have also seen that (62) is equivalent to F(X, p) = 0, which allowed us 
to solve for A as a function of p and finally obtain the equation Q{p) = 0, where Q : TL —¥ l\. By Lemma 27, 
1Z = kcrCf'(O) © A4, and Q'ifi) : Ai —> If is a linear homeomorphism. Hence for p G 1Z we may write p = £ + 77, 
where £ G ker^'(O) and 77 G M.. By the implicit function theorem, we can solve for 77 in terms of £, and this 
equation defines a smooth manifold in a neighbourhood of with infinite dimension and codimension, since the 
spaces kerC/'(0) and M are infinite-dimensional by Lemma 27. □ 



8 Conclusions and open problems 

In this paper, we considered the planar bilaplacian with a smooth, radially symmetric and compactly supported 
potential 9 and described the set of perturbations of the potential in the space 1Z = L 2 ([0, ri]; H 1 / 2 (5' 1 ), r dr) 



39 



for which an embedded eigenvalue persists. We expect that the space 1Z can be replaced by the Sobolev space 
H 1 / 2 {B ri (0)) of H 1 ' 2 -functions of two variables that have support in the ball B ri (0). 

One restriction of our work is that we consider only potentials with compact support: The reason is that we 
were forced to work with different function spaces of solutions for r small and for r large. For r small, we have 
some freedom in choosing the space, as any space of the form X = H 1+a x H a x H 1 x L 2 , with < a < 1, 
ensures that an exponential dichotomy exists. For r large, due to the structure of the equations, there is no 
such freedom, the regularity on the first two components has to be same as the regularity of the last two. So 
it is unclear whether there exists an exponential dichotomy when the support of p is not compact. It would be 
interesting to see whether our hypothesis that p has compact support could be replaced by an appropriate decay 
condition on p. 

For the original potential 9, we see no obstacles in removing the condition that 9 has compact support. It should 
be possible to replace this condition by the long range condition |#'(r)| < C(l + r)~ 1- ^ for some /3 > 0. It should 
also be possible to remove the condition that 9 is radially symmetric, although considerably more work will be 
needed without this condition. 

We believe that the methods put forward in this paper can be used to study other operators. In particular, the 
exponential-dichotomy results established in [18] arc for systems of reaction-diffusion equations, so we believe 
that the only obstacle for extending our results to selfadjoint systems are the presence of nonsmooth potentials. 
For other operators, it might not be possible to modify the function spaces involved to prove the existence of 
exponential dichotomies. These are difficult problems that have to be studied in future work. 
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